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Resumen 


En este trabajo práctico integrador correspondiente al curso de Mecánica del Vuelo, se 
llevó a cabo un estudio exhaustivo y análisis detallado de la aeronave asignada al grupo, el 
Short 360. El enfoque se centró en investigar aspectos clave como la performance de vuelo, su 
estabilidad estática y estabilidad dinámica. Se procedió a realizar comparaciones meticulosas 
entre los valores obtenidos en el análisis y aquellos proporcionados por el fabricante, además 
de efectuar una validación rigurosa en concordancia con las normativas pertinentes. 


Palabras clave: Performance, estabilidad estática, estabilidad mecánica, Short 360. 


Abstract 


In this integrated practical work for the course Flight Mechanics, an in-depth study and 
analysis were conducted on the group-assigned aircraft, the Short 360. The focus encompas- 
sed exploring critical aspects including flight performance, static stability, and dynamic stabi- 
lity. Comprehensive comparisons were drawn between the analysis results and manufacturer- 
provided data, accompanied by meticulous validation in alignment with relevant regulations. 


Keywords: Performance, static stability, dynamic stability, Short 360. 
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1. Introducción 


En el presente documento se abordará la elaboración colaborativa del proyecto integrador gru- 
pal para la asignatura de especialización en Mecánica del Vuelo. En los contenidos de esta carpeta 
se explorarán y analizarán los aspectos relacionados con la Performance, Estabilidad Estática y 
Estabilidad Dinámica de la aeronave asignada al grupo, el Short 360. 


1.1. Descripción de la aeronave 


El Shorts 360, también conocido como SD3-60, es una aeronave de transporte regional desa- 
rrollada por la compañía británica Short Brothers plc. Este proyecto surgió en respuesta al creci- 
miento del transporte aéreo regional en la década de 1970 y la demanda de aerolíneas europeas de 
aeronaves de este tipo. En 1980, se presentó al público el desarrollo del Short 360. El 1 de junio 
de 1981, el primer prototipo realizó su vuelo inaugural. Después de tres meses, el 3 de septiembre 
de 1981, la aeronave recibió su certificado tipo 2019). El 19 de agosto de 1982, el primer 
avión de producción completó su primer vuelo y en noviembre de 1982, comenzó a operar con 
Suburban Airlines (Corner) [2023). En la Figura[1]se muestra el prototipo con matrícula G-ROOM. 


Figura 1: Short 360 


Se construyeron un total de 165 unidades hasta 1991 y se muestra una visualización en 3 vistas 


de la aeronave en la Figura [2] obtenida del manual de vuelo de la aeronave (Short Brothers plc. 
1985). 
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Figura 2: Visualización en 3 vistas de la aeronave asignada 
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Se presentan algunas características importantes de la aeronave en la siguiente Tabla[I] 


Tabla 1: Características de la Aeronave 


Características generales 


Designación Simbología Valor Unidad 
Envergadura b 22.80 m 
Superficie alar S 42.73 m? 
Altura hiso: 7.21 m 
Longitud Lio: 21.58 m 
Cuerda raíz C> 2.27 m 
Cuerda puntera C; 1.92 m 
Peso en vacío W+;acio 7870.00 kg 
Peso máximo Woax 12000.00 kg 
Tripulación Norew 3 - 
Capacidad Noariencer 36 - 
Perfil raíz NACA 634418 
Perfil puntera NACA 634414 
Perfil E.H. NACA 0012 
Perfil E. V. NACA 0016 
Motor Pratt € Whitney Canada PT6A - 65R 1062 kW 
Hélice Hartzell Model HC-B5MP-3C/¡MI0876ASK 5 palas 


App Joe 
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2. Performance 


2.1. Introducción 


La performance de una aeronave engloba una serie de atributos que definen el comportamiento 
y las capacidades de un avión en relación con su propósito de fabricación. Estos atributos esta- 
blecen límites en magnitudes críticas, como velocidades máximas y mínimas, tasas máximas de 
ascenso y otros parámetros esenciales. 

La adecuada determinación y comunicación de estos aspectos resulta de suma importancia para 
garantizar una Operación segura y eficiente de la aeronave. Al entender y tener pleno conocimiento 
de las capacidades inherentes de la aeronave, la tripulación puede tomar decisiones informadas du- 
rante el vuelo y adaptarse a diferentes situaciones. Esto es especialmente crucial en situaciones de 
emergencia o condiciones climáticas adversas, donde el conocimiento preciso de las capacidades 
de la aeronave puede marcar la diferencia entre un aterrizaje seguro y una situación peligrosa. 

Como principio fundamental en la ingeniería aeronáutica, estas capacidades deben estar en 
línea con las regulaciones y normativas establecidas por las autoridades de aeronavegabilidad. 
Estas normativas no solo garantizan la seguridad de las operaciones aéreas, sino que también 
buscan mantener un nivel de calidad y rendimiento uniformes en la industria de la aviación. 

En este sentido, este trabajo se enfocará en explorar en profundidad los diversos aspectos que 
componen la performance de una aeronave, examinando cómo estas características influyen en 
la operación y en la experiencia de vuelo. Además, se buscará verificar el cumplimiento de las 
normativas de aeronavegabilidad vigentes que aplican al tipo de aeronave que se está estudiando. 


2.2. Vuelo estacionario, recto y nivelado 


2.2.1. Introducción 


En esta sección, procederemos a llevar a cabo un análisis exhaustivo de la performance du- 
rante el vuelo nivelado de la aeronave asignada al grupo, el Short 360. El enfoque principal de 
esta sección consiste en determinar, mediante una metodología detallada que expondremos pos- 
teriormente, la denominada "envolvente de vuelo” de la aeronave. Esta curva distintiva define de 
manera precisa el rango dentro del cual la aeronave es capaz de mantener un vuelo nivelado y 
recto. Además de presentar la envolvente, se exhibirán las gráficas correspondientes a la potencia 
requerida y la potencia disponible en diversas altitudes, tomando en consideración la variación de 
la velocidad. 


2.2.2. Metodología 


Para el desarrollo del cálculo se siguió la metodología presentada por la cátedra, la cual se 


desarrolla en plenitud en la bibliografía recomendada (Roskam y Lan| 1997). 


2.2.2.1 Hipótesis Simplificativas 


Las hipótesis utilizadas para la metodología propuesta, se detallan a continuación: 
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= Se considera una aeronave en vuelo nivelado estacionario donde los parámetros de velocidad 
y altitud permanecen constantes. 


= Ángulos de incidencia y trayectoria recta pequeños y = 0 vuelo recto. 


= La aeronave se encuentra en condición trimada, es decir se considera la ausencia de momen- 
tos y fuerzas netos 


= Se asume que todas las fuerzas actúan en el centro de gravedad de la aeronave y que la masa 
de la aeronave es constante todo el tiempo. 


= Se considera la ausencia de velocidades angulares 


2.2.2.2 Datos necesarios 


Los datos necesarios se presentan en la Tabla [2] los mismos fueron obtenidos por medio de 
relevamientos y del manual de vuelo de la aeronave ” Flight Manual for the SD3-60 ” (Short 


Brothers plc.||1985). 


Tabla 2: Datos necesarios 


Designación Simbología Valor Unidad 
Superficie alar Su 42.17 m? 
Peso de la aeronave W 94019 N 
Coef. de sustentación máximo Cz 1.7763 - 


Velocidad de pérdida 


Para determinar el valor de C,,,,., fue necesario obtener la velocidad de pérdida. Para este 


propósito, se utilizó el gráfico representado en la Figura [3| el cual fue extraído del manual del 
fabricante 11983). 

Dado que la velocidad de pérdida obtenida del gráfico en la Figura[3|corresponde a la velocidad 
calibrada (CAS), se procede a calcular la velocidad verdadera (TAS). Es importante señalar que se 
asume un factor f. igual a uno, ya que no se tienen en cuenta los efectos de compresibilidad. 


0,5 
EAS = f.-CAS =TAS l db ) 
PMSsL 
TAS=Y E 46|"| 
Ss 
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Figura 3: Gráfica Peso vs. Velocidad de pérdida 


Potencia Requerida 


Siguiendo las hipótesis previamente establecidas y considerando las fuerzas que actúan en una 
condición de vuelo recto y nivelado, podemos derivar las ecuaciones generales de movimiento 
para esta situación: 


T=D 


L=W 


Dado que disponemos de una curva polar, detallada en el Apéndice en la Tabla que se 
puede expresar de la siguiente manera: 


Cp = Cp, + AC, + BC? (1) 


Donde los valores son los siguientes: 


Cp, = 0,0261 
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A = 0,0090 


B = 0,0357 


Podemos derivar la expresión para calcular la Potencia Requerida P,.,, que se define como se 
muestra en la Expresión |2] 


1 1 
Preg = TV = DV = ¿PS vc) = ¿PS v? (Co, + ACL:+ BC;) Q) 


Potencia Disponible 


El manual del fabricante de la aeronave y el manual del fabricante del motor (Ver Ref. 
Brothers plc.! [19851 Pratt 8% Whitney|2007)) proporciona información sobre las curvas de potencia 


del motor a diferentes altitudes (presentadas como torque del motor). Estos datos se recopilaron 
numéricamente para permitir el cálculo de las potencias disponibles (P,,,). El detalle de las curvas 
de potencia se encuentra en el Apéndice [6.2 

A partir de esta información, podemos derivar la expresión para la potencia en el eje del motor: 


2 
Peje =2-T E - 1,3558 - 1300) [Watts] 3) 


en donde el torque se define de la siguiente manera: 


T = -0,0897h + 3625 [lb-ft] (4) 
Considerando el rendimiento de la hélice, podemos calcular la Potencia Disponible P,,, de la 
siguiente manera: 


Pos = Mn * Peje (5) 


Eficiencia de la Hélice 


La eficiencia de la hélice, la cual varía con la velocidad, se calcula numéricamente siguiendo 
el enfoque propuesto por Solies. La curva de eficiencia resultante se presenta en el Apéndice, en 
la Expresión[81] junto con una gráfica en la Figura[65|ilustrativa de los resultados del método y el 
ajuste numérico propuesto. 

La expresión obtenida es la siguiente: 


Tm ==2,7837x 10?-V*+1,0935x 10*.V?*-1,6159x10*.V?+1,1873x107?.V+4,3272x 10 
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Velocidades Mínimas y Máximas 


Las velocidades mínimas y máximas se obtienen cuando la potencia requerida es igual a la po- 
tencia disponible. Matemáticamente, esto corresponde a igualar las expresiones P]y[] Al hacerlo, 
se llega a la siguiente ecuación: 


2B -Ww? Cm-:5:-VW-p 
Ss V-p 2 


Esta ecuación es de cuarto orden y sus raíces representan las velocidades mínimas y máximas 
para diferentes altitudes de vuelo: 


Prep = A: V-W+ = Po, 


Co Pos A. w.y-p yr 2B WE 


> ET (6) 


2.2.3. Resultados 
2.2.3.1 Curvas de Potencia disponible Vs Potencia requerida 


Una vez ajustadas las curvas de potencia en el eje y de la eficiencia de la hélice, es posible 
obtener las curvas de potencia disponible empleando la Expresión [5] Las curvas que describen la 
potencia requerida se obtienen a partir de la ecuación dada en la Expresión|2] Estas curvas, tanto la 
de potencia requerida como la de potencia disponible, han sido trazadas para diferentes altitudes y 
se presentan conjuntamente en la Figura/4] 


2.2.3.2 Envolvente de Vuelo 


Mediante un procedimiento numérico, se calcularon los puntos de intersección de las curvas 
de potencia requerida y potencia disponible de la aeronave, es decir, los puntos en los que ambas 
potencias son iguales y se obtienen las velocidades máximas y mínimas. Estos puntos se han 
recopilado y se presentan gráficamente como la envolvente de vuelo de la aeronave, que se muestra 
en la Figura[5] 

La ecuación a resolver se encuentra dada por la Expresión[7] 


Pieg(v, h)=- Pay(h) =0 (7) 


Para llevar a cabo la resolución numérica, se utilizó la librería scipy y su módulo optimize. 
Se aplicó una versión modificada del método de Powell diseñada para calcular raíces múltiples 
(consultar |scipy.optimize.root). 

Además de la representación gráfica de la envolvente de vuelo en la Figura [5] se trazaron las 
velocidades de pérdida para diversas altitudes. Estos valores se obtuvieron utilizando la ecuación 
que hace uso del coeficiente máximo de sustentación Cz 


max * 


2 
V, = ia (8) 
S Cha? 
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Figura 4: Potencia requerida vs Potencia disponible. 
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Mon 
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Figura 5: Envolvente de Vuelo de la Aeronave 
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2.2.4. Discusión 


En esta sección, se han obtenido representaciones gráficas de la potencia requerida y la poten- 
cia disponible a diferentes altitudes (ver Figura (4). Estas curvas han permitido definir una envol- 
vente de vuelo que delimita una región específica en la cual la aeronave puede mantener un vuelo 
recto y nivelado (ver Figura|5). 

El techo absoluto calculado corresponde a una altitud de 6269 metros, en comparación con 
el techo absoluto de la aeronave especificado como 7010 metros (23000 pies) en (doc8643.com! 
2023). Se observa una discrepancia significativa de 741 metros, que equivale a un error relativo 
porcentual del 10.57 %. Estas diferencias podrían atribuirse a las simplificaciones de la teoría 
utilizada y posiblemente a imprecisiones en el cálculo de la curva polar de resistencia. 

En cuanto al rango de velocidades obtenido, los resultados son coherentes con las velocidades 


óptimas de vuelo presentadas en el manual de la aeronave (Short Brothers plc. (1985). En ese 


sentido, se considera que los resultados obtenidos son satisfactorios en términos generales. 
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2.3. Vuelo en trepada estacionario 


2.3.1. Introducción 


En esta sección, se llevará a cabo el análisis de la etapa de vuelo en trepada de la aeronave 
asignada a nuestro grupo, el Short 360. La fase de vuelo en trepada es fundamental para que la ae- 
ronave alcance su altitud de crucero y pueda cumplir su función de transporte de carga o pasajeros. 
Se determinarán los gráficos de la tasa de ascenso y el gradiente de ascenso a diversas altitudes 
de vuelo para la aeronave. Asimismo, se calculará el valor del techo de servicio de la aeronave y 
se superpondrán los valores de velocidad para la tasa de ascenso máxima en diferentes altitudes 
dentro de la envolvente de vuelo definida en la sección anterior. También se determinará el techo 
de servicio con un motor inoperativo. Por último, se verificará el cumplimiento de las regulaciones 
de aeronavegabilidad correspondientes con respecto a los requisitos de vuelo en trepada. 


2.3.2. Metodología 


Para el desarrollo y la determinación de lo descrito anteriormente se seguirá la metodología 
propuesta por la cátedra la cual se encuentra desarrollada en plenitud en la bibliografía (Roskam y 


1997 [Ojhal (1995). Se parte de las siguientes hipótesis: 


= Vuelo en trepada siguiendo una trayectoria de vuelo aproximadamente recta y a ángulos (y) 
de ascenso relativamente bajos (menores a 157). 


= Vuelo en trepada con velocidad de la aeronave constante, es decir vuelo en ascenso estacio- 
nario. 


2.3.3. Resultados 
2.3.3.1 Tasa de ascenso RC y gradiente de ascenso CGR 


Se aplicó el procedimiento iterativo recomendado en la bibliografía y por la 
cátedra para determinar los parámetros de ascenso específicos de la aeronave asignada. Como 
resultado de este proceso, se generaron las curvas correspondientes a la tasa de ascenso (RC) y al 
gradiente de ascenso (CGR) de la aeronave. 

El valor del gradiente de ascenso, (CGR), representa la cantidad de altitud ganada por unidad de 
distancia horizontal recorrida por la aeronave. Posteriormente, este valor se empleó para calcular 
las distintas tasas de ascenso en función de la altitud. Las representaciones gráficas de estas curvas 
se encuentran en las Figuras[6] y [7] respectivamente. 


Proceso iterativo seguido 


El proceso iterativo seguido es el mencionado anteriormente y se detalla a continuación: 


1. Sean la altitud h, velocidad V, peso W, superficie alar S, curva polar Cp(C,) y potencia 
disponible P,,, prescritos. 
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2. Se empieza por proponer un gradiente de ascenso inicial sin(y) = 0. 


3. Se obtiene el coeficiente de sustentación correspondiente de: 


_ Weos(y) 
== 10V2S 


4. De la curva polar se obtiene el coeficiente de resistencia correspondiente a ese Cz. 


5. Se obtiene la potencia requerida para la resistencia aerodinámica en la trepada de la forma: 


1 
Pa=DWVe ¿VS Cp 
6. Se obtiene entonces un nuevo gradiente de ascenso de la forma: 


(Pay => P+eq) 


sin(y) = WV 


Se reemplaza luego este nuevo valor en el paso 3 y se repite hasta la convergencia. 


7. Una vez obtenido el gradiente de ascenso se obtiene la tasa de ascenso con la ecuación: 


RC = V sin(y) 
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Figura 6: Gradiente de ascenso para distintas altitudes 
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Figura 7: Tasa de ascenso para distintas altitudes 
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2.3.3.2 Obtención del techo de servicio de la aeronave 


Utilizando los valores de tasa de ascenso obtenidos y representados en la Figura[7| empleamos 
Python para calcular los valores de RC, en diferentes altitudes de vuelo. El proceso iterativo 
proporcionó soluciones discretas para CGR y RC, y estos valores se almacenaron en vectores O 
listas. A continuación, se utilizaron funciones integradas de Python, como max(), para determinar 


los RC, , COrrespondientes a cada conjunto de datos de altitud (consulte Built-in Functions). 


Con los valores de RC, , en Mano, se generó el gráfico que se muestra en la Figura][8] 


mM Valores de RCmas obtenidos 
—— Ajuste lineal, h = —16.1611 - RCimax + 6287.9570 | 
*  RCmaw para techo de servicio (FAR25) | 


0 50 100 150 200 250 300 350 400 
ROmax - ([m/min] 


Figura 8: Tasas de ascenso máximas 


Luego, se realizó un ajuste lineal a los puntos obtenidos, que se muestra en la Expresión [S] 
Siguiendo las normativas aplicables a la aeronave en análisis (FAR 25), se determinó que el techo 
de servicio se alcanza a una tasa de ascenso de RC ax = 200 [ft/min] = 60,96 [m/min], y su valor 
se describe en la Expresión |10] Del mismo modo, para una tasa de ascenso de RC ja, = O [m/min], 
se obtiene el techo absoluto, cuyo valor se presenta en la Expresión [1 1] 


h.=—16,1611 + RCar + 02879370 (9) 
Mserr =3302,7776 [Mm] =3303 [ma] (10) 
Rabs =06237,9570 [m] = 6288 [mi] (1) 
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2.3.3.3 Obtención del techo de servicio OEI 


Se realiza de forma similar a lo que se detalló en la sección anterior pero en este caso se tuvo 
en cuenta la mitad de la potencia disponible. Obtenidos los valores de RC, se generó el gráfico 
que se muestra en la Figura[9] 


mM Valores de RCmas obtenidos 
—— Ajuste lineal, h = —27.8160 + RCimax + 3394.2036 | 
*  RCmax para techo de servicio (FAR25) 


0 20 40 60 80 100 120 140 
RCmazw - [m/min] 


Figura 9: Tasas de ascenso máximas - OEI 


Luego, se realizó un ajuste lineal a los puntos obtenidos, que se muestra en la Expresión 
Siguiendo las normativas aplicables a la aeronave en análisis (FAR 25), se determinó que el techo 
de servicio se alcanza a una tasa de ascenso de RC ax = 200 [ft/min] = 60,96 [m/min], y su valor 
se describe en la Expresión [13] Del mismo modo, para una tasa de ascenso de RC ax = 0 [m/min], 
se Obtiene el techo absoluto, cuyo valor se presenta en la Expresión [14] 


h =-27,8160 - RCnax + 3394,2036 (12) 
hservo,, = 1698,5392 [m] = 1699 [m] (13) 
habso:y = 3394,2036 [m] = 3394 [m] (14) 
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2.3.4. Discusión 


En esta sección, se han calculado y representado gráficamente las tasas de ascenso y los gra- 
dientes de ascenso para diversas altitudes de vuelo. Las Figuras [7] y [6] muestran estas representa- 
ciones respectivamente. 

Luego, utilizando un método que implica ajustes lineales junto con extrapolaciones simples, 
se determinaron el techo de servicio y el techo absoluto de la aeronave. Para este análisis, se con- 
sideró una tasa de ascenso máxima de referencia de 200 ft/min, que es un estándar para aeronaves 
certificadas bajo la norma FAR 25. Los resultados obtenidos fueron los siguientes: el techo de ser- 
vicio se estimó en 5303 m, mientras que el techo absoluto, que corresponde a una tasa de ascenso 
nula, fue de 6288 m. 

Adicionalmente, se realizaron los mismos cálculos para la condición de un motor inoperativo 
(OEI, por sus siglas en inglés), considerando la mitad de la potencia disponible. Los resultados 
obtenidos en esta situación fueron los siguientes: el techo de servicio con un motor inoperativo se 
estimó en 1699 m, y el techo absoluto con un motor inoperativo se calculó en 3394 m. 

En relación a las regulaciones de aeronavegabilidad vigentes, se ha clasificado esta aeronave en 
función de sus características dentro de los parámetros establecidos por el FAR 25, considerando 
su naturaleza como aeronave de transporte. Específicamente en lo que respecta a su desempeño de 
ascenso, se establece que el gradiente de ascenso (CGR) para una aeronave bimotora en fase de 
trepada en configuración limpia (es decir, con los flaps retraídos y el tren de aterrizaje en posición 
recogida y asegurada) no debe ser inferior a 0.012, según lo estipulado en la referencia 
2014). Los resultados obtenidos se ajustan satisfactoriamente a estos requisitos normativos, lo que 
indica que la aeronave cumple con los estándares establecidos. 

Finalmente, se incorporaron estos resultados en la envolvente de vuelo de la aeronave, como 
se muestra en la Figura 
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Figura 10: Superposición de los resultados obtenidos en la envolvente de vuelo de la aeronave 
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2.4. Vuelo en maniobra 


2.4.1. Introducción 


En esta sección, se procede a detallar el desempeño de la aeronave Short 360 en maniobras 
de giro coordinado. Se entiende por maniobra una operación que implica al menos una rotación y 
cambio de dirección. Estas rotaciones se logran mediante la deflexión de las superficies de control 
de la aeronave, como los alerones (para el rolido), el estabilizador horizontal (para el cabeceo) y 
el timón de dirección (para la guiñada). 

La condición de vuelo asociada a un giro coordinado se define como aquella en la cual la 
aeronave tiene un ángulo de inclinación ¿ con respecto a la línea del horizonte, y su trayectoria de 
vuelo sigue una órbita circular con un radio de giro R, en el plano horizontal. 

A continuación, se detalla el procedimiento para calcular diversos parámetros relevantes en la 
performance del giro coordinado, como el radio de giro, la velocidad de giro y el factor de carga. 


2.4.2. Metodología 


Para el desarrollo del cálculo se siguió la metodología presentada por la cátedra, la cual se 


desarrolla en plenitud en la bibliografía recomendada (Roskam y Lan||1997| 2023d). 


2.4.2.1 Factor de carga máximo en giro coordinado 


El factor de carga máximo en esta condición de vuelo se expresará en función de la velocidad, 
que se encuentra en el rango de velocidades desde la velocidad de pérdida hasta la velocidad 
máxima para vuelo recto y nivelado. Además, este factor de carga también variará según la altitud 
en la que se encuentre la aeronave. 

Es importante señalar que el factor de carga se encuentra limitado por diversos aspectos, que 
incluyen consideraciones aerodinámicas, la potencia disponible y la resistencia estructural. Por 
lo tanto, se procederá a determinar el factor de carga máximo en función de cada una de estas 
limitaciones. 

El factor de carga máximo limitado por consideraciones aerodinámicas, denotado como Pax (Cima)? 
se calcula según la Expresión|15] que representa la relación entre el coeficiente de sustentación má- 
ximo de la aeronave C,,,,, y el coeficiente de sustentación de la aeronave en vuelo recto y nivelado 
a la velocidad seleccionada Ch see: 


_ Comas 
Amax,Cimax el 


(15) 
Cr 


El factor de carga máximo limitado por la potencia disponible se calcula mediante la Expresión 


CLPa 
MmaxPa = HT (16) 
Ciel 
Donde Ch», se determina a partir de la curva polar, previamente hallando el coeficiente de 
resistencia máximo Cp,,,, con la Expresión [17] 
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C Pas (17) 
Dmax»Pav FET 
5pVv3S 
Finalmente, el factor de carga máximo na se obtiene como el valor mínimo entre los factores 
de carga limitados por la potencia disponible, la resistencia estructural y la aerodinámica: 


Nmax = MIA (Mmax.c Lmax? Amax,P ay > act) 


en donde el factor de carga estructural Max.s1 tiene un valor de Mmax.sir = 2,/38, el mismo se 
encuentra dado en el manual del fabricante y también en el certificado tipo de la aeronave (Short 


Brothers ple][1985)JEASA] (2019). 


2.4.2.2 Radio de giro 


Una vez encontrado el factor de carga máximo, se procede a calcular el radio de giro R, para 
la altura a la que se encuentra la aeronave en función de su velocidad empleando la siguiente 
expresión: 


2.4.2.3 Velocidad de giro 


Obtenido el radio de giro, es posible calcular la velocidad de giro o velocidad angular de 
guiñada y mediante la siguiente relación cinemática: 
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2.4.3. Resultados 


Utilizando la metodología detallada previamente, se ha calculado el valor del factor de carga 
máximo (Na) en el rango de velocidades comprendido entre la velocidad de pérdida y la velocidad 
máxima de la aeronave para diversas altitudes. Es importante destacar que para velocidades fuera 
de este rango, el factor de carga será inferior a la unidad, lo que impide el cálculo del radio de giro 
y la velocidad de giro. 

Una vez obtenidos los valores del factor de carga, se procedió a calcular y graficar el radio de 
giro en el mencionado rango de velocidades para distintas altitudes. Los resultados se presentan 
en la Figura [11] Asimismo, se realizó el cálculo y la representación gráfica de las velocidades de 
giro en el mismo rango y a diferentes altitudes, lo que dio lugar a la Figura[12] 

Se generaron gráficos del factor de carga correspondiente a la condición de vuelo en crucero 
en el rango de velocidades, destacando cada una de sus limitaciones: aerodinámicas, de potencia 
y de resistencia estructural. Este gráfico se presenta en la Figura 

Finalmente, en las Figuras [14] y se muestran las curvas que representan el radio de giro 
mínimo y la velocidad de giro máxima en función de la altitud de vuelo, respectivamente. 
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h=0 [m] 
h = 1000 [m] [ 
h = 2000 [m] 
2000 h = 3000 [m] 
h = 4000 [m] 
h = 5000 [m] 
[ 


h = 6000 [m] 


— 1500 


1000 


R; - [m 


500 


50 60 70 80 90 100 110 
V - [m/s] 


Figura 11: Radios de giro para distintas altitudes 
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Figura 12: Velocidades de giro para distintas altitudes 
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Figura 13: Limitaciones de factor de carga para la condición en vuelo crucero 
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Figura 14: Radio de giro mínimo para la aeronave Short 360 
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Figura 15: Velocidades de giro máximas para la aeronave Short 360 


28 de [126] 


Desimone - Alvarez Granda 


Mecánica del Vuelo: UNC ES pedia 
Trabajo Integrador Grupal - Short 360 


2.4.4. Discusión 


En esta sección, se llevaron a cabo los cálculos de los factores de carga correspondientes a las 
diferentes limitaciones de la aeronave. Es importante destacar que para este tipo de avión, el factor 
de carga estructural está limitado a n = 2,758. Por lo tanto, los resultados obtenidos indican que 
las limitaciones por factor de carga provienen exclusivamente de la aerodinámica y la potencia 
disponible en la planta propulsora. 

Una vez calculados los factores de carga para diversas altitudes, se generaron las gráficas pre- 
sentadas en las Figuras [11]y [12] Estas representaciones gráficas muestran la variación de los radios 
y las velocidades de giro en función de la velocidad de vuelo, manteniendo una altitud constan- 
te. Además, las curvas se componen claramente de dos ramas: una relacionada con la limitación 
aerodinámica y otra asociada a la potencia disponible. 

Los gráficos de las Figuras [14] y proporcionaron información sobre la evolución del radio 
de giro mínimo (R,,,,) y la velocidad de giro máxima (Ymax) en función de la altitud de vuelo. 
Los resultados observados fueron coherentes, ya que a medida que la altitud aumentó, se registró 
un incremento en el valor de R,,,, y una disminución en la velocidad de giro máxima Ya. Este 
comportamiento es característico en aeronaves, ya que a altitudes más elevadas, las limitaciones 
aerodinámicas en los factores de carga se vuelven más restrictivas. En consecuencia, para realizar 
giros coordinados a estas altitudes, la aeronave requerirá un radio de giro más amplio y, por ende, 
deberá ejecutar esta maniobra a una velocidad menor. 
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2.5. Alcance y autonomía 


2.5.1. Introducción 


En esta sección, se abordará el análisis de la capacidad de alcance y autonomía de la aeronave 
asignada, el Short 360. Para ello, se considerará una variación en el peso de la aeronave, que estará 
directamente relacionada con las fluctuaciones en el consumo de combustible. 

El objetivo principal de esta sección es determinar el consumo de combustible de la aeronave 
en diferentes condiciones, incluyendo la fase de trepada, el comportamiento durante el vuelo en 
crucero a una altitud y velocidad constantes, y finalmente, la creación de un diagrama de carga 
paga vs. alcance que represente los resultados obtenidos. 


2.5.2. Metodología 


Para llevar a cabo el cálculo, hemos seguido la metodología establecida por la cátedra, la cual 


se encuentra detallada en su totalidad en la bibliografía recomendada (Roskam y Lan||1997). En 


esta fase, abordaremos nuevamente el equilibrio en las ecuaciones de movimiento, pero ahora 
incorporaremos la variación de peso debida al consumo de combustible. 


2.5.2.1 Datos 


Los datos utilizados para la obtención del alcance y de la autonomía son los que se muestran 


en la Tablal3] 


Tabla 3: Datos 


Datos iniciales 


Parámetro Símbolo Valor Unidad 
Peso máximo al despegue Wurow 11657 kg 
Peso en vacío operacional Wacio 7480 kg 
Carga paga máxima WePne: 3102 Kg 
Carga de combustible máxima Weomimáx 1742 kg 
Reserva de 45 minutos [FR Was min 281 kg 
Velocidad crucero Vencer 100 m/s 
Altitud de crucero Mercero 3050 m 
Densidad a la altitud de crucero  Perucero 0.9044 kg/m? 
Consumo específico Cp 1.17E-06  N/w/s 


"El consumo específico en crucero fue estimado en el Apéndice|6.4 
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2.5.3. Resultados 


2.5.3.1 Combustible consumido en la trepada 


Para determinar el consumo de combustible durante la fase de trepada, se empleó un método 
iterativo recomendado por la cátedra, el cual se encuentra descrito en mayor detalle en la biblio- 
grafía (Saldía]/2023b) 

En las condiciones de trepada, se estableció un ascenso desde el nivel del mar hasta la altitud 
de crucero de 3050 [m] (un valor dentro de la envolvente de vuelo). Se consideraron tanto trepadas 
con la aeronave en su peso máximo al despegue (Wurow) como el escenario en el que la aeronave 
trepa sin carga paga y con la máxima cantidad de combustible disponible. Para determinar el 
consumo de combustible durante la trepada, se resolvió la siguiente integral: 


hoy, 
f.climb 
Wrclimb = f dh 
hi 


dh ! 
en donde RC = — es la tasa de ascenso y W  ctimp es el flujo de combustible. 


Se obtuvo una altura intermedia entre las alturas inicial al nivel del mar y la de crucero que son 
hi; =0my horucero = 3050 m respectivamente, esto es: 


Se realizó lo mismo para obtener un peso intermedio, partiendo desde el Wyrow a partir de la 
siguiente expresión: 


We ctimb 

2 

Luego se obtuvo la tasa de ascenso RC para la condición del peso W;_» y altura h,-2 y velocidad 
de trepada v;rep = 77 [m/s] constante. 

Se obtiene entonces el flujo de combustible de la siguiente expresión, hasta que se logró la 
convergencia: 


W¡-2 = Wurow — 


Wi. 
O E E 
f.climb ( 1 2) RC 49 


Se obtuvieron los siguientes resultados: 


= Para el caso de trepada con la aeronave trepando con Wyrow se obtuvo: 


a = 68,51 [Kg] 


= Para el caso de trepada con la aeronave trepando sin carga paga y con máximo combustible, 
se realizó lo mismo, pero en vez del peso Wyrow se utilizó W acio + We, y se obtuvo: 


Wireradas = 43,50 [Kg] 
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2.5.3.2 Pesos al inicio y final de crucero en cada condición 


Los pesos al inicio y final del crucero en cada condición son los siguientes: 


= En el caso de la condición de carga paga máxima mas combustible, el peso al inicio 
del crucero se calcula de la siguiente manera. En primer lugar, se determina el peso del 
combustible para esta condición siguiendo: 


Weóniba = Wurow = Wer. > Weacio = 1075 [kg] 


Por lo tanto el peso al inicio y al final del crucero para ésta condición se obtienen de: 


W; => Meacio + Web: + Wicoriba a Wena = 11588 [Kg] 


W, = W; a Wesmba + Was min = 10862 [Kg] 
= Para la condición de máximo combustible mas carga paga, se obtiene primero el peso de 


la carga paga: 
Wer, = Wurow = Wieemmnás E Woacio = 2435 [kg] 


Por lo tanto el peso al inicio y al final del crucero para ésta condición se obtienen de: 
W = Wacio + Weom,max + Wer, => Wirepada = 11588 [kg] 
W, => Wacio uE Wer, + Was min = 10195 [kg] 


= Para la condición de máximo combustible sin carga paga, el peso al inicio y al final del 
crucero es: 
W; = Wiacio + Wenas = Weradar = 9179 [kg ] 


W) = Woiacia + Was min ="7761 [Kg] 


2.5.3.3 Confección del diagrama de carga paga vs alcance 


Para determinarlo, se realizaron cálculos de alcance en diversas condiciones de crucero, man- 
teniendo una altitud y velocidad constantes (h — V, a una altitud de 3050 [m] y una velocidad de 
100 [m/s]). En todos los casos, se consideró una reserva de combustible equivalente a 43 minu- 
tos de vuelo adicionales (Considerando que la aeronave realiza sus vuelos de forma IFR 
2000). 

El alcance y la autonomía se obtuvieron mediante la resolución de las siguientes integrales 
respectivas: 


Wa Wa 
Mp Cz 1p Ch 
R=-| BT gy im E=- | == ies] 
le Cc) Cp W wm Cp V Cp W 8 


Los resultados obtenidos se listan a continuación: 
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= Para la condición de carga paga máxima más combustible 


>| R;,-y=385,83 [km] y Ep-y = 1,07 [h] 


= Para la condición de máximo combustible más carga paga 


>| Ry =1353,00 [km] y E-y =3,76 [h] 


= Para la condición de máximo combustible sin carga paga 


>| R;,y=1562,57 [km] y E,-y = 4,34 [h] 


Finalmente el diagrama de carga paga Vs. Alcance se muestra en la Figura|16| 
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Figura 16: Diagrama de Carga paga vs Alcance para la Aeronave Short 360 
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2.5.4. Discusión 


En esta sección, se llevaron a cabo cálculos para determinar la capacidad de alcance y autono- 
mía de la aeronave asignada. Inicialmente, se calculó el consumo de combustible durante la fase de 
trepada utilizando el método iterativo recomendado por la cátedra, en dos condiciones diferentes. 
La primera considera la aeronave en su peso máximo al despegue, lo que resultó en un valor de 
Wiepada, = 68,51 Kg, y la segunda considera la capacidad máxima de combustible sin carga paga, 
obteniendo Wrepadaz = 43,50 Kg. 

Luego, se calcularon los pesos al inicio y al final del crucero para cada una de estas condicio- 
nes: carga máxima de combustible más carga paga, máxima capacidad de carga paga más combus- 
tible y máxima capacidad de combustible sin carga paga. Estos cálculos se realizaron manteniendo 
la aeronave a una altitud y velocidad constantes (3050 m y 100 m/s). 

Posteriormente, se generó el gráfico de carga paga versus alcance, que se muestra en la Figu- 
ra Se consideró una reserva de combustible para 45 minutos de vuelo, considerando que la 
aeronave realiza su vuelo utilizando las reglas de vuelo instrumentales (IFR) (Ver 2000). 

Además, se calcularon las autonomías para cada una de estas condiciones, obteniendo los si- 
guientes resultados: 1,07 horas de vuelo para la condición de carga paga máxima más combustible, 
3,76 horas de vuelo para la condición de máxima capacidad de combustible más carga paga y 4,34 
horas de vuelo para la máxima capacidad de combustible sin carga paga. 

Se compararon los resultados obtenidos con los datos proporcionados en la publicación de 
Jane's All the World's Aircraft (ver (Taylor|[1981)). Según esta fuente, la aeronave tiene un alcance 
de 425 [km] al volar a una altitud de crucero de 3050 [m] (10000 [f£]) a una velocidad de 200 
nudos (103 [m/s]) con carga paga máxima mas combustible. En nuestro cálculo dicha condición de 
crucero arrojo un resultado de 386 [km], los resultados obtenidos son próximos a los que comenta 
la publicación, si bien hay diferencias podemos considerar que son dentro de todo aceptables. 
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2.6. Despegue y aterrizaje 


2.6.1. Introducción 


En esta sección, analizaremos el rendimiento de despegue y aterrizaje de la aeronave asignada 
a nuestro grupo, el Short 360. El objetivo principal de esta parte del estudio es determinar la 
longitud total necesaria para el despegue y aterrizaje de la aeronave bajo condiciones específicas 
de peso, altitud y temperatura en la pista de despegue. 

Este análisis es esencial para garantizar un desempeño seguro y eficiente de la aeronave. Ade- 
más, es crucial cumplir con las regulaciones de seguridad aérea vigentes, ya que esta información 
influye en las limitaciones de operación de la aeronave en distintos aeropuertos. 


2.6.2. Metodología 


La condición de despegue y aterrizaje constituye un caso particular dentro del ámbito de las 
ecuaciones de movimiento, ya que implica un vuelo con aceleración. En este contexto, se consi- 
derará la configuración de flaps deflectados y se utilizará la curva polar de baja velocidad, además 
de asumir que no hay viento. 

El cálculo se basa en la metodología proporcionada por la cátedra y se encuentra detallado de 


manera exhaustiva en la bibliografía recomendada (Roskam y Lan|/1997). 


Para el caso que abordaremos a continuación, se han tomado las condiciones correspondientes 
al Aeropuerto Córdoba / Ing. Aer. Ambrosio Taravella (CBA - SACO - COR), el cual se encuentra 
a una elevación de 489 metros sobre el nivel del mar y se considera un ángulo de inclinación de 
pista nulo (4 = 0). 


2.6.2.1 Distancia de despegue 


Para obtener la distancia total de despegue se fueron calculando las distancias de cada etapa 
que comprende el mismo, cómo se presenta en la Figura[17] 


Distancia de despegue 


SrTo 


Distancia en el aire 


SA 


Distancia de rodaje 
Sa 


Nariz en pista Trepada 


Transición 


Rotación 


Trepada_.-- 


-” 


Robetáculo 


Figura 17: Fases de despegue 


Dónde el cálculo de cada etapa se detalla a continuación: 


Desimone - Álvarez Granda 35 de 


Trabajo Integrador Grupal - Short 360 


= Distancia de Nariz en pista S yar: 
Para obtener ésta distancia se resolvió la siguiente integral: 


V 
* V 
Suoe= | dd V 
0 


Donde el valor del coeficiente de sustentación teniendo en cuenta el efecto suelo, se obtuvo 


de la siguiente expresión obtenida de 1940): 


1 
Ch, = FA he 


= Distancia de rotación S p: 
Se obtiene de resolver la siguiente ecuación, considerando un tiempo de rotación tg = 3 seg 
ya que es un avión de transporte: 


1 
Sr= > (Ve + Vzor) tr 


= Distancia de transición S Tp: 
Para obtener ésta distancia se calcularon primeramente los parámetros que intervienen, los 
mismos son: 


Coeficiente de sustentación adicional: 


1 [(ViorY Vs Y 
AC; = — —1 —- 38 


_ 2W/S 
"E" pgAC; 
Obtención del ángulo de trayectoria: 
. T-D 
sin(y.) = y 


Finalmente: 
STR = Rrp sin(y.) 


= Distancia de trepada S c;: 
Para obtener esta distancia, primeramente se obtuvo la altura al final de la transición, donde 
Msercen = 35 ftes la altura libre de obstáculos obtenida de la norma FAR 25: 


hrr = Rrkr — Rrx cos (y4) 
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siendo la distancia en la trepada: 


an > MrR 
tan (y.) 


ScL = 


En el caso en donde la altura de transición hy resulte mayor a la altura libre de obstáculos 
Mscreen» la distancia correspondiente a la trepada S ¿, no será tenida en cuenta. 


La distancia total requerida para el despegue es la suma de las distancias correspondientes a 
cada etapa, es decir: 
Sro=SNGrR+FSRESTRES CL 
2.6.2.2 Distancia de aterrizaje 


Análogamente que con la distancia de despegue, se calculan las distancias de cada etapa que 
comprende el aterrizaje, las mismas se muestran en la Figura 


Distancia de aterrizaje, Sy 


Distancia aterrizaje en el suelo Distancia aterrizaje en el aire 


SLa LA 


Nariz en el suelo Rotación Transición Descenso 


Obstáculo 


Figura 18: Fases de aterrizaje 


Dónde el cálculo de cada etapa se detalla a continuación: 


= Distancia en fase de descenso S ¡pes Se calculan los parámetros que intervienen en su ob- 
tención, esto es: 


Angulo promedio de la aproximación: 


Altura al inicio del flare: 
Hare = KR lar (1 cos Ya) 


Por lo tanto la distancia de esta etapa es: 


iva ll Flare 


S LDES = tany 
A 
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= Distancia en fase de transición Í ¿rr: 
Se obtiene primeramente el radio del flare, esto es: 


2 
Vitare 
R flare A E 
g (A siare > 1) 
Y la distancia recorrida será: 
SLTR = RflareYA 


= Distancia con todas las ruedas en contacto con el suelo S ¿nar: 
Esta distancia se obtiene de resolver la siguiente integral: 


0 V 
a f e 
dd Vrp 5 (T Ñ D, — Mgbrake (w Ni Le) + 0,08 W (teams Ñ s)) 


Finalmente la distancia total de aterrizaje se obtiene de la suma de las fases mencionadas 
anteriormente, esto es: 


SD =3S DES + SLTR + S LNGR 
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2.6.3. Resultados 


Los datos utilizados para la determinación de las distancias de despegue y aterrizaje se presen- 
tan en la Tabla l4P?] 


Tabla 4: Datos para cálculo de distancias de despegue y aterrizaje. 


Datos Simbología Valor Unidad 
Superficie alar S ref 42.17 mé 
Peso de despegue MTOW 12000.00 kg 
Peso al aterrizaje MLW 11840.00 kg 
Coeficiente de sustentación máximo al despegue Cua 2.4333 - 
Coeficiente de sustentación máximo al aterrizaje ¿A 2.7065 - 
Flaps para despegue ÓTo 15 2 
Flaps para aterrizaje ÓLD 30 d 
Coeficiente de rozamiento con la pista He 0.0250 - 
Coeficiente de rozamiento debido al frenado Hbrake 0.8000 z 
Coeficiente de sustentación con efecto suelo Ch, 0.3823 - 
Ángulo de inclinación de la pista 10) 0.00 d 
Factor de carga en el flare ilar 1.08 - 
Polar para despegue Cp = 0,0511 -— 0,0187C, + 0,0360C% 
Polar para aterrizaje Cp = 0,0590 — 0,0216C, + 0,0360C? 

Velocidades para el despegue 
Velocidad de pérdida V, 44.24 m/s 
Velocidad de rotación Va = 1,10V, 48.66 m/s 
Velocidad de despegue Vior = 1,15V, 50.88 m/s 
Velocidad al final de la fase de trepada V, = 1,20V, 53.09 m/s 
Velocidades para el aterrizaje 

Velocidad de pérdida V, 41.67 m/s 
Velocidad de aproximación V, = 1,30V, 54.17 m/s 
Velocidad en el flare Vs = 0,95 Va 51.46 m/s 
Velocidad de contacto con la pista Vip = 1,15V, 47.92 m/s 


2.6.3.1 Distancia de despegue 


A continuación se presentan los resultados obtenidos para las distancias de despegue con AEO 
(Todos los motores operativos) y OEI (Un motor inoperativo) en las Tablas [5] y [6] 


2Las velocidades de pérdida junto a los C,,,. para el despegue y aterrizaje fueron determinados en el Apéndice|6.5 
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Tabla 5: Resultados obtenidos para la distancia de despegue para AEO 


Distancia de despegue AEO 


Parámetro Símbolo Valor [m] 
Distancia con nariz en pista  SyGR 325 
Distancia de rotación Sr 149 
Distancia de transición STR 584 
Distancia de trepada S cr 0.00 
Distancia total de despegue S ro 1058 


Tabla 6: Resultados obtenidos para la distancia de despegue para OEI 


Distancia de despegue OEI 


Parámetro Símbolo Valor [m] 
Distancia con nariz en pista  SyGr 738 
Distancia de rotación Sr 149 
Distancia de transición STR 134 
Distancia de trepada ScL 192 
Distancia total de despegue STo 1213 


Finalmente, la distancia total para el despegue resulta el valor mayor de aceleración despegue 


con OEI y 1.15 aceleración despegue con AEO (FAR 25.113 - Transport Category Airplanes: 
1998). Los valores obtenidos para ambos casos se presentan en las Expresiones[18]y [19] 


>| Srowo =1,15- 1058 [m] = 1217 [m] (18) 
>| Sroo, = 1213 [m] (19) 


Tomando como se dijo anteriormente el máximo de ambos valores, la distancia de despegue 
total resulta lo dado en la Expresión [20] 


> STo = 1217 [1] 


Para la resolución de las integrales se empleó el método de integración de Simpson con un 
dv = 0,010, el mismo se aplicó utilizando Python a través de la librería Scipy con su módulo 
Integrate. 


(20) 


2.6.3.2 Distancia de aterrizaje 


A continuación en la Tabla[7]se presentan los resultados obtenidos en el cálculo de la distancia 
total de aterrizaje. 
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Tabla 7: Resultados obtenidos para la distancia de aterrizaje 


Distancia de aterrizaje 


Parámetro Símbolo Valor [m] 
Distancia de descenso S LDES 151 
Distancia de transición (flare) S LTR 438 
Distancia con nariz en pista S LNGR 175 
Distancia total de aterrizaje S 1D 764 


La distancia total de aterrizaje resulta lo dado por la Expresión |21] 


Para la resolución de las integrales se empleó el método de integración de Simpson con un 
dv = 0,010, el mismo se aplicó utilizando Python a través de la librería Scipy con su módulo 
Integrate. 


2.6.3.3 Velocidad de decisión 


Se determinó la velocidad de decisión (V¡) mediante el criterio de pista balanceada (o equili- 


brada) (Saldía) [2023c). Los resultados que se obtuvieron se muestran en la Figura[19] 


—— Distancia Aceleración - Parada 


—— Distancia Aceleración - Despegue ' 


O Vip crítica 


0 b) 10 15 20 25 30 39 40 45 50 
Ver = [m/s] 


Figura 19: Criterio de pista balanceada 
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Resolviendo numéricamente mediante Python se determinó la intersección de ambas curvas y 
se obtuvo el valor de la velocidad de decisión o crítica; el resultado de la misma se presenta en la 


Expresión|22] 
> V, = 46,28 [m/s] (Q2) 


2.6.3.4 Datos provistos por el fabricante 


Para calcular la distancia de despegue, recurrimos al manual del fabricante (Short Brothers plc. 
1985). Siguiendo el algoritmo proporcionado en dicho manual y considerando las características 


del aeropuerto seleccionado como referencia, el Aeropuerto Internacional de Córdoba Ingeniero 
Aeronáutico Ambrosio Taravella (FAA: CBA; IATA: COR; OACI: SACO), se detallan a continua- 
ción las condiciones clave para hacer uso del algoritmo del manual del fabricante: 


= Temperatura ambiente: 11,8215"C (aproximadamente 12'C) 
= Peso de despegue: MTOW = 12000 kg 


= Altitud de presión del aeropuerto: 3298 [ft] (aproximadamente 3300 [ft]). Se determinó 
utilizando una referencia de ONH según la atmósfera estándar, que es de 955,86 Hpa, y 
considerando que el aeropuerto se encuentra a una elevación de 1604 [ft] (489 m) con 


respecto al nivel medio del mar (consultar Figura[22] (ANAC|/2023)). 


Con estos parámetros mencionados anteriormente, el resultado proporcionado por el fabricante 
para la distancia de despegue es de aproximadamente 3500 [ft] (aproximadamente 1067 [m]), 
como se muestra en la Figura|20] 

Para el caso de la distancia de aterrizaje se mantuvo el aeropuerto y las condiciones, solo 
cambió el peso de la aeronave el cual para esta condición se adoptó el peso máximo posible 
para el aterrizaje MEW = 11840 kg, asumiendo condiciones de viento en reposo y declinación 
de pista nula. El resultado proporcionado por el fabricante para la distancia de aterrizaje es de 
aproximadamente 2100 [ft] (aproximadamente 640 [m]), como se muestra en la Figura/21] 
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Figura 20: Distancia de despegue del manual del fabricante 


Desimone - Álvarez Granda 43 de [126] 


Mecánica del Vuelo: E r” pa AA 
Trabajo Integrador Grupal - Short 360 de Aeronáutica 


ALTERNATE LANDING FIELD LENGTH 


<< 3838.22 1332. 0.00.b20(ÁA.xM.0 


o 
A 7] 001/sauLam 
A 0001/1334 


0 ap 
RUNWAY ORADIEN? 


2 1 
DOWN 


10 20 30 
HEAD 


0 
TAIL 
REPORTED WIND 
SPEED-KNOTS 


10 


3 


. 
O" FREUD 


GAUSITELANIDSNS 
> TAC 


LANDING FIELD LENGTH 
(UNFACTORED) 


20 21 22 23 24 25 26 27 


AIRCRAFT WEIGHT 
LB/1000 


O Lo +++ +] 
17 18 19 


canasonos 


10 


PRESSURE 
FEET/1000 


5 


ALTITUDE- 


0 
AIRFIELD 


vYO'9S9UA YES" va 


Figura 21: Distancia de aterrizaje del manual del fabricante 
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Figura 22: Detalle del Aeropuerto Internacional de Córdoba Ingeniero Aeronáutico Ambrosio Ta- 
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2.6.4. Discusión 


En esta sección, se han calculado las distancias requeridas tanto para el despegue como para el 
aterrizaje de la aeronave bajo estudio, obteniendo valores de 1217 metros y 764 metros, respecti- 
vamente. En términos generales, consideramos que estos valores son adecuados para la aeronave 
analizada, ya que se sitúan dentro del rango esperado para una aeronave de transporte regional. 

Es común que una aeronave turboprop regional, como la que estamos evaluando, necesite en- 
tre 1500 y 2000 metros de pista para operar eficazmente. Estas cifras le otorgan la flexibilidad de 
Operar en una amplia variedad de aeropuertos, desde pequeños hasta medianos. A modo de ejem- 
plo, la aeronave podría operar sin dificultades en aeropuertos como el de Córdoba, que cuenta con 
pistas de 3200 metros en la pista 01/19 y 2280 metros en la pista 05/23, el Aeroparque Metropo- 
litano Jorge Newbery, que dispone de una pista de 2100 metros en la pista 13/31, o el Aeropuerto 
Internacional Teniente Luis Candelaria de San Carlos de Bariloche, que cuenta con una pista de 
2200 metros en la pista 29/11, por mencionar algunos ejemplos de aeropuertos nacionales. 

Es relevante destacar que el fabricante también proporciona información para calcular las dis- 
tancias de despegue y aterrizaje en las condiciones de máximo peso de despegue (MTOW) y 
máximo peso de aterrizaje (MLW), que son 1067 metros y 640 metros, respectivamente, según se 
detalla en (ver Figuras [20] y [21). Es importante tener en cuenta que 
los valores calculados pueden diferir de los proporcionados por el fabricante debido a las posi- 
bles variaciones en las condiciones exactas utilizadas en los cálculos. Además, es plausible que 
existan diferencias sustanciales en la estimación de la curva polar de baja velocidad para ambas 
condiciones y esto haya inducido cierto error en los cálculos. 
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3. Estabilidad Estática 


3.1. Introducción 


El propósito fundamental de llevar a cabo un análisis de estabilidad estática radica en la minu- 
ciosa evaluación de las fuerzas y momentos involucrados. Esto se realiza con el fin de garantizar 
que, en presencia de una perturbación, estas magnitudes tiendan a realinear la aeronave hacia su 
posición de equilibrio original. Esto implica que la suma de todas las fuerzas y momentos que 
actúan sobre la aeronave debe mantener la condición de equilibrio. Este equilibrio se traduce en 
la conservación del vuelo trimado, que es el estado en el cual la aeronave mantiene su actitud y 
trayectoria deseada sin necesidad de correcciones constantes por parte del piloto. 


3.1.1. Objetivo principal del análisis de estabilidad estática 


= Mantener la condición de vuelo trimado (vuelo en equilibrio alrededor de algún eje). 


= Lograr un equilibrio de fuerzas (lo cual ya se estudiaba en el análisis de performance) y 
momentos (lo cual se agrega en este análisis). 


= Analizar las fuerzas y momentos para determinar si tienden a mantener la condición de 
equilibrio, es decir, si restauran el equilibrio después de una perturbación. 


3.2. Estabilidad Estática Longitudinal 


3.2.1. Introducción 


En esta sección, se llevará a cabo el análisis de la Estabilidad Estática Longitudinal de la aero- 
nave Short 360 en estudio. El objetivo principal es determinar los márgenes del centro de gravedad, 
encontrando las posiciones del punto neutro con y sin potencia, con mando fijo y con mando li- 
bre, tanto para vuelo nivelado como para vuelo en maniobra. También se buscará identificar las 
posiciones más adelantadas del centro de gravedad permisibles para lograr el vuelo trimado. 

Además, se analizarán otros parámetros relevantes, como la curva de sustentación para el vuelo 
trimado y la curva de ángulo de elevador requerido para el trimado de la aeronave en función de la 
velocidad de vuelo. 


3.2.1.1 Hipótesis 


= Se asume un movimiento simétrico en el plano de simetría de la aeronave, donde los despla- 
zamientos están restringidos a este plano. Se considerará únicamente la rotación alrededor 
del eje Y, (cabeceo) (Ver Figura[23). 


= El empuje se ajusta para equilibrar la fuerza de resistencia mediante la posición del acelera- 
dor del motor. 
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= La sustentación se ajusta para equilibrar el peso, utilizando un coeficiente de sustentación 
CL que varía según la velocidad y la altitud de vuelo. 


= El momento de cabeceo M., alrededor del centro de gravedad, causado por el ala, el fuse- 
laje, la planta motora y el empenaje, debe ser equilibrado por el momento generado por el 
elevador. 


Plano de simetría 


Figura 23: Aeronave, su plano de simetría y su sistema de ejes cuerpo 


3.2.1.2 Objetivos 


= Determinar los puntos neutros para diferentes condiciones de vuelo. Estos puntos son de 
particular interés, ya que representan las posiciones límite del centro de gravedad para las 
cuales se obtiene una derivada C,,,, negativa (indicativa de estabilidad) o positiva (indicativa 
de inestabilidad). 


= Calcular los márgenes estáticos para cada condición de vuelo. 


Los datos relevantes que se utilizarán en esta sección se presentan en la Tabla[8| 


3Para ver la estimación de algunos de los parámetros mas relevantes, consultar el Apéndice|6| 
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Tabla 8: Datos necesarios - Estabilidad Estática Longitudinal 


Parámetro Símbolo Valor Unidad 
Posición del centro aerodinámico ala-fuselaje boe 0.0611 Adim. 
Relación de presiones dinámicas Mi 0.9500 Adim. 
Pendiente de sustentación E.H. a; 4.3523 1/rad 
Pendiente de sustentación ala - fuselaje Dub 6.2290 1/rad 
Superficie alar S 42.7300 m? 
Superficie de E.H. Sa 9.8300 m? 
Superficie del elevador Se IIDID m? 
Gradiente de Downwash 0€e/0a 0.3640 Adim. 


0.0597 rad 
6.8687 l/rad 
0.2650 Adim. 


Downwash a incidencia nula € 

Pendiente de la configuración completa a 

C.G. de la aeronave h 

Distancia entre C.A. ala - fuselaje y E.H. Í, 10.9850 m 

Distancia entre C.A. empenaje horizontal y C.G. del avión L, 10.5377 m 

Cuerda Aerodinámica Media Cc 1.9900 m 

Volumen de cola Vi 1.2182 Adim. 
Vi 


Volumen de cola 1.2699 Adim. 
Diámetro de la Hélice D, 2.8200 m 
Ángulo de pala B 30.0000 ñ 
Área de la hélice 5 6.2458 m? 
Distancia desde el plano de la hélice hasta 1/4 CAM la 1.8613 m 
Coeficiente de cabeceo libre del ala-fuselaje Crca..; -0.2106 Adim. 
Ángulo de calaje del E.H , -4.5000 2 
Pendiente de sustentación del perfil E.H (NACA 0012) Ch, 5.7869 1/rad 
Pendiente de sustentación mando libre a! 6.5640 1 /rad 
Cuerda aerodinámica media del E.H. Us 1.5800 m 
Masa de la aeronave m 9587.2739 kg 


3.2.2. Margen estático para mando fijo con y sin potencia 


En la condición de Mando Fijo, las superficies de control se mantienen en una posición fi- 
ja y controlable, lo que significa que el parámetro de control, la deflexión del elevador, no está 
influenciado por fuerzas aerodinámicas, sino que se ajusta de manera que permanezca en una po- 
sición constante. El margen estático se define como la diferencia entre las posiciones del centro de 


gravedad (4) y el punto neutro (/h,,) 2005; /Etkin y Reid] [1996). 
K, =h,-h 


El punto neutro con mando fijo /h, se obtiene cuando C,,,, = 0, lo que indica una estabilidad 
neutra. Su cálculo se presenta en la Expresión|23] 
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(23) 


10/04 


Oe€ 1 0C my 
a 0% 


An = An, O m1 


Para la condición sin potencia, se ignora la contribución del momento de la fuerza normal 


debido al sistema propulsivo a la estabilidad longitudinal, representado por 
calculó de acuerdo con la siguiente expresión (Saldía|[2023a): 
Cm Ss » Xp Oe | 0C»,, 
ra SS da 


Mp > Al 
. Este término se 
100% 


00) 


0€, 
En donde el gradiente de upwash -—— = 0,30 fue estimado utilizando el Reporte NACA WRL25 
a alas elípticas (consultar Ref. on Ribner) [1944)) a través de la Figura [24] Además, el término 


7e. =0,27 1/rad, que refleja la variación del coeficiente C,,, con el ángulo de incidencia sobre el 

ca de la hélice «,, se estimó utilizando el Reporte NACA 819 (consultar Ref. (Ribner) [1945)) 
para una relación de avance de J = 1,64 y un ángulo de pala de 6 = 30”. Dado que la aeronave 
cuenta con una hélice de 5 palas, se tomó un promedio entre los gráficos correspondientes a hélices 
Hamilton Standard de 6 y 4 palas, como se muestra en la Figura [25] 


0 0 0 0 79 O O O O O O O O A 
A 
A MAA A 
MES O 
MEMO 
ME O 
000000 


£.0 


EFE 
Rua 
AS 


Diatence behiná root quarter-ohord point, in root chorda 


Figure 5,- Value of 1 - Er on longitudinal axie of elliptic wing 
for aspeot ratios 6, 9, and 12, 


Figura 24: Estimación del gradiente de Upwash 
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FIGURE 8. —Slde-force derivative for singlerotating Hamilron Standard propeller 3185-48 wlth <pínner. Slr blades; 0.182, 


Figura 25: Obtención de la derivativa para fuerza normal 
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Los puntos neutros para mando fijo, tanto sin como con potencia, resultaron: 


> A, =0,57531 > | Ah, =0,3635 


Asp Mcp 


Por otro lado los valores de los margenes estáticos para mando fijo resultaron: 


>| K,, =0,3101 >| K,, =0,2985 


3.2.3. Margen estático para mando libre con y sin potencia 


En la condición de mando libre, las superficies de control están libres, y la posición que adopten 
dependerá de las fuerzas aerodinámicas. Esto implica que el coeficiente de momento de charnela 
asume un valor de (Cy,) = 0, y que el ángulo del elevador resulta ser el llamado ángulo de flotación 
de = Oe yr 

Se define al margen estático con mando libre como la diferencia entre las posiciones del centro 
de gravedad (h) y el punto neutro con mando libre (»/,), expresado como: 


K;, = (lr, —h) 
z . . / Cno Ch, 
El cálculo del punto neutro con mando libre se obtiene al plantear C,,, = Cin, — e =0 
He, 
(estabilidad neutra). La expresión para calcular h/, se presenta en la ecuación de la Expresión|24] 
Ch,, Ch, 1 a Cy 
Ml, =h, + hy, - — (a.Vy + Ch, An z 24 
n a'Ch,, a (a H Ese e Ch, ( ) 


o ., 00 Ñ ; ; 
en donde la variación A = Cz,, del coeficiente de sustentación de la configuración completa por 
e 
la deflexión del elevador se calculó como: 


S S 
Ch, = ma = MATE = 0,5975 


en donde el valor de la efectividad del elevador 7 = 0,59 se obtuvo con la Figura [26] para una 
Se 
relación de superficies — = 0,38. 


Í 
La estimación de los coeficientes de charnela con efectos tridimensionales se realizó mediante 
el procedimiento recomendado por la cátedra y el mismo se encuentra detallado en el Reporte de 


la NACA WR L - 663 (consultar Ref. 1943)). Se utilizaron las Figuras y Los 


coeficientes de charnela tridimensionales resultaron luego de aplicar el método lo siguiente: 


Ch, =-0,0010 [1/DEG] Cy, =-0,0045 [1/DEG] 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 
S./S; - [Adim.] 


Figura 26: Obtención del parámetro de efectividad del elevador 7 
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Figura 27: Reporte de la NACA WR L - 663 


Los puntos neutros para mando libre, tanto sin como con potencia, resultaron: 
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» Pr, = 0,4721 » Pr, = 0,4603 


Por otro lado los valores de los margenes estáticos para mando libre resultaron: 


>| K;, =0,2071 >| K;, =0,1953 


3.2.4. Puntos neutros de maniobra para mando fijo y mando libre 


En el caso de vuelo recto nivelado estacionario, se requiere un ángulo del elevador ó, y una 
fuerza en el mando P para mantener la condición de equilibrio de momento longitudinal. En el 
caso de realizar una maniobra, se necesitarán cambios Aó, y AP para mantener el equilibrio en 
esta condición de vuelo, ya que aparece una velocidad angular de cabeceo (q) que modifica los 
coeficientes aerodinámicos de sustentación y momento (AC; AC,,,). Dado que asumimos el vuelo 
en maniobra como estacionario (sin aceleración angular), se tiene AC,, = O (para mantener el 
equilibrio). 


AS, _ CwC1, (24-C,) 


n= 2uA 
La posición del centro de gravedad para la cual el ángulo del elevador por la gravedad es nulo 
(2) = 0, se denomina punto neutro de maniobra con mando fijo y se define como se muestra 


en la ecuación de la Expresión[26| 


(A — hn) (25) 


h, =h Cn 26 
m — " 2u4=C;, ( ) 


Donde el factor de densidad relativa uu se calcula como: 


_ 2m 


 = — = 249,3359 
psc 
Las derivadas C,, Ci,» y Ch.,> que indican la dependencia de los coeficientes aerodinámicos 
de sustentación y momento, y el coeficiente de momento de charnela, con la velocidad angular de 


cabeceo (q), se calculan como: 


Ci, = (C,,), = 29:0,Vy = 10,6518 [1/rad] 
l 
Cy % (Cry), = 21:41 Vu ==56,4045 [1 /rad] 


/ 
Cu, =2Chn, = = -2,1123 [1/rad] 


eq a = 
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Estos cambios en los coeficientes aerodinámicos son debidos al incremento del ángulo de 
incidencia sobre el empenaje horizontal debido a la velocidad angular (q) en la maniobra. Por otro 
lado, el punto neutro de maniobra con mando libre está dado por la ecuación de la Expresión 


A [€ Ch, 
a A ) 


=_——— 27 
Ch, Ch. Cr, 2u - Ch, ( ) 


Donde la variación del coeficiente de momento de charnela por el ángulo de ataque de la 
configuración completa (a) se obtuvo como: 


= Ch, 1 - 7) = 0,0364 [1/rad] 
di OQ 


Finalmente, los puntos neutros de maniobra, para mando fijo y mando libre, sin y con potencia, 
resultaron respectivamente: 


Ch 


04 


>| hn, = 0,6881 >|, =0,6093 


>| hn, =0,6765 » a = 0,5977 


y los márgenes estáticos: 


>| K,, =0,4231 >| K;, =0,3443 


>| Km, =0,4115 » daa =0,3327 


3.2.5. Diagrama de puntos neutros y C.G. 


Con los resultados obtenidos en las secciones anteriores, se ha creado el gráfico de la Figura|28] 
el cual muestra las posiciones de los puntos neutros para las condiciones de mando fijo y mando 
libre, considerando y sin considerar los efectos de potencia. En cuanto al vuelo en maniobra, los 
puntos neutros se han representado de manera similar a lo detallado anteriormente, como se puede 
apreciar en la Figura Finalmente, en la Figura se presentan los puntos neutros para las 
condiciones de mando fijo y mando libre tanto en vuelo nivelado como en vuelo en maniobra. 
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—— h,, M.F. Sin Potencia 
—— h,, M.F. Con Potencia 
—— |, M.L. Sin Potencia 
— |, M.L. Con Potencia 
—— h posición C.G. 


0.0 0.2 0.4 0.6 0.8 1.0 
x/c - [Adim.] 


Figura 28: Posición de los puntos neutros de mando fijo y mando libre 


—— him M.F. Sin Potencia 
—— hn MF. Con Potencia 
—— h, M-L. Sin Potencia 
— hi, M.L. Con Potencia 
—— h posición C.G. 

0.0 0.2 0.4 0.6 0.8 1.0 


x/c - [Adim.] 


Figura 29: Posición de los puntos neutros de mando fijo y mando libre para vuelo en maniobra 
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Figura 30: Posición de los puntos neutros para vuelo nivelado y vuelo en maniobra 
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3.2.6. Curva de sustentación en vuelo trimado 


Para entrar en la condición de vuelo trimado, es necesario que se cumpla: 


Ci =0 


Al definir un coeficiente de sustentación de trimado para el vuelo deseado, se llega a un sistema 
de dos ecuaciones con dos incógnitas que relacionan dicho coeficiente con el ángulo de elevador 
necesario para esta condición. Luego de solucionado es posible obtener la ecuación para el coefi- 
ciente de sustentación en la condición de trimado, el mismo se encuentra dado por la Expresión 


(28) 
m5 

En donde el coeficiente de momento total de cabeceo en relación con el ángulo de ataque se 
obtiene mediante la siguiente expresión: 


Cm, = Ch, (h 5 Mryp) 0 a + sl 
' : da da 
Cuando se desplaza el elevador, se produce una modificación en el ángulo de ataque de todo el 
empenaje horizontal, lo que a su vez afecta la sustentación y el momento alrededor del centro de 
gravedad de la aeronave. 
La capacidad del elevador para modificar el coeficiente de momento se conoce como potencia 
del elevador C,,,,, y se define como: 


5 S 
Cmo = Mide | Va — ls — hos) 


En la Figura [31]se muestra gráficamente la pendiente de sustentación para vuelo trimado sin 
potencia. 
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Figura 31: Curva de sustentación en vuelo trimado 


3.2.7. Curva de ángulo del elevador necesario para mantener el vuelo trimado 
El ángulo del elevador necesario para vuelo trimado se calcula a través de la Expresión [29] 


do. 
== Ó. + trim C : 29 
0 dC. Lerim ( ) 


El primer término, correspondiente a la sustentación nula se obtiene de: 


Ctrim 


e Cu Cr. 


eo A 
Donde el determinante A se define como: 


o) 


A= Er Cm = Ca Ets 


Por último, la derivada de la posición del elevador con respecto al coeficiente de sustentación 
para vuelo trimado se obtiene de: 


Ae 
AC Lin o 


Se realiza la curva para la condición crucero para el valor de h = 0,265 dado en la Tabla [8] 
teniendo y no teniendo en cuenta los efectos de potencia, esto mismo se encuentra representado en 


Ch 
“(h-h, 
e ) 


Ctrim 
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la Figura[32] Luego en la Figura[33|se realiza el mismo gráfico solo que en función de la velocidad 
equivalente. 


— der. - Sin Potencia 
—— Se, - Con Potencia 


0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 
Ch - [Adim.] 


Figura 32: Curva de ángulo del elevador necesario para h = 0,265 
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—— der, - Con Potencia 
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Figura 33: Curva de ángulo del elevador necesario para h = 0,265 con velocidad equivalente 


3.2.8. Curva de Ángulo de Elevador Requerido por g 


Se ha obtenido la curva que representa el ángulo de elevador requerido por g en función de la 
velocidad equivalente, utilizando la condición de vuelo crucero con potencia. Para este propósito, 
se empleó la Expresión[25|(Repetida mas abajo). La representación gráfica de esta curva se muestra 


en la Figura[34] 


A9.  CwC.,(2u-C,,) e 
n=t 24 A (M1 = hm) 


Desimone - Álvarez Granda 61 de [126] 


Trabajo Integrador Grupal - Short 360 


60 70 80 90 100 110 120 
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Figura 34: Curva del Ángulo de Elevador Requerido por g en función de la Velocidad Equivalente 


3.2.9. Posición mas adelantada del C.G. 


Se determinó la posición más adelantada del centro de gravedad teniendo en cuenta la deflexión 
máxima del elevador obtenida del certificado tipo (EASA|[2019) y considerando el C,, calculado 


max 


en Performance. Este cálculo se basa en la pendiente de la recta que se obtiene al considerar los 
puntos de ó., con C, =0 y Ós,,, = —12* con C, = C,,, Con esta pendiente, podemos despejar la 


max 


posición más adelantada del centro de masas utilizando la fórmula de la Expresión [30] 


max 


dó. A 
dC; a 


» ta = 0,1186 (31) 


Luego sumando este resultado al gráfico dado en la Figura[30]|se obtiene el gráfico de la Figura 


h=h,=- 


(30) 


de esto último se obtuvo: 
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hn, M.F. Sin Potencia 

h, M.F. Con Potencia 
h/, M.L. Sin Potencia 

h?, M.L. Con Potencia 
hm M.F. Sin Potencia 
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Figura 35: Puntos neutros junto a la posición de h mas adelantada. 


3.2.10. Condición de despegue 


Finalmente teniendo en cuenta la condición de despegue de la aeronave y aplicando nueva- 
mente la metodología anteriormente detallada, se determinaron nuevamente las posiciones mas 
adelantadas y atrasadas posibles del centro de gravedad, los resultados obtenidos se representan 
en la Tabla [9] y en la Figura [36] 


Tabla 9: Resultados obtenidos para la condición de despegue 


Resultado Símbolo Valor 
Punto Neutro de Mando Fijo Sin Potencia An, 0,5739 
Punto Neutro de Mando Fijo Con Potencia ep 0,5656 
Punto Neutro de Mando Libre Sin Potencia A, 0.4482 
Punto Neutro de Mando Libre Con Potencia Po. 0.4397 
Punto Neutro de Maniobra con Mando Fijo Sin Potencia Asp 0.6693 
Punto Neutro de Maniobra con Mando Fijo Con Potencia Amo 0,6609 
Punto Neutro de Maniobra con Mando Libre Sin Potencia Pts 0.5798 
Punto Neutro de Maniobra con Mando Libre Con Potencia Pep 0.3715 
Posición mas adelantada Madelmax 0.2106 
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Figura 36: Puntos neutros junto a la posición de h mas adelantada en condición de despegue. 


3.2.11. Discusión 


En esta sección, se ha llevado a cabo un minucioso análisis de la Estabilidad Estática Lon- 
gitudinal de la aeronave asignada al grupo, el Short 360. Los resultados obtenidos han permitido 
visualizar cómo el factor de la potencia introduce inestabilidad como también las condiciones de 
Mando Libre, que se revelaron como las más críticas en términos de estabilidad, tanto en vuelo 
crucero como en condiciones de despegue. Como conclusión, se confirma que, tal como se espera- 
ba para una aeronave de este tipo, el Short 360 muestra un comportamiento estáticamente estable 
en todas las condiciones de vuelo analizadas en este apartado. 

Por otro lado, las curvas de ángulo de elevador en función del coeficiente de sustentación en 
vuelo trimado presentan un comportamiento dentro de los parámetros normales. Se observa un 
desplazamiento hacia abajo de estas curvas, así como una disminución de la pendiente cuando se 
quita el factor correspondiente a la potencia. Del mismo modo, las curvas de ángulo de elevador 
en función de la velocidad equivalente de vuelo (EAS) también muestran un desplazamiento hacia 
abajo cuando se quita el factor correspondiente a la potencia. 
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3.3. Estabilidad Estática Direccional 


3.3.1. Introducción 


En esta sección, se llevará a cabo un análisis de la estabilidad direccional de la aeronave asig- 
nada al grupo, el Short 360. El enfoque se centrará en el aspecto estático de la estabilidad con 
respecto a perturbaciones laterales contenidas en el plano horizontal. Se supone que la respuesta 
de la aeronave es lo suficientemente lenta como para desestimar los efectos de la velocidad de 
guiñada y la velocidad de rolido. 

En la presente sección se empleará la metodología presentada en la bibliografía recomendada 


por la cátedra 2004; U.S. Department of Defense] [1977). Los datos requeridos para esta 


sección se encuentran dados en la Tabla[10] 


Tabla 10: Datos necesarios - Estabilidad Estática Direccional 


Parámetro Símbolo Valor Unidad 
Pendiente de sustentación E. V. a, 4.2288 1/rad 
Ángulo de flecha al cuarto de cuerda Aoja 0.3015 5 
Ángulo de diedro r 3.0000 S 
Alargamiento A 12.1657  Adim. 
Distancia entre el CG y el CA o Xa 0.2039 Adim. 
Relación de presiones dinámicas Nv 1.0000 Adim. 
Volumen de cola del E. V. v, 0.0833 Adim. 
Superficie alar $5 42.7300 m? 
Superfice lateral del fuselaje S Bs 44.5453 m? 
Superficie del E.V. e 7.7300 m? 
Envergadura b 22.8000 m 
Longitud del fuselaje ls 21.2000 m 
Longitud del fuselaje / 4 15/4 5.3000 m 
Longitud del fuselaje - 3/4 lp-3/4 — 15.9000 m 
Altura del fuselaje a 15/4 hi 2.8092 m 
Altura del fuselaje a l5- 34 h> 2.1366 m 
Reynolds del Fuselaje Res 1.13E+08 — Adim. 
C.G. desde la nariz del fuselaje pa 8.5301 m 
Deflexión del rudder máxima a 221 R 
Efectividad del timón T 0.62 Adim. 
Velocidad de crucero v. 100.0000 m/s 
Velocidad de despegue Viof 50.8800 m/s 
Velocidad de aproximación Vo 54.1700 m/s 
Densidad de crucero De 0.9044 kg/m? 
Desnidad en la aproximación Pa 1.1685 kg/m? 
Coef. Resistencia Crucero Cb. 0.0316 Adim. 
Coef. Resistencia Aterrizaje Cb, 0.1294 Adim. 
Distancia del motor al plano de simetría d 3.1000 m 
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3.3.2. Criterio de estabilidad 


Una aeronave es direccionalmente estable si esta es capaz de re-alinearse con la resultante del 
viento, siempre y cuando sea perturbada desde un vuelo nivelado. 


Ng >0 
Cr > 0 
En donde el coeficiente de momento de guiñada está dado por: 
N 
Ca ==: 
asb 
0C, 
Crg = 
B 06 


Se asume que C, es la suma de las contribuciones individuales causadas por el fuselaje, las 
alas y las superficies de la cola, cómo así también las contribuciones debida a la planta propulsora. 


3.3.2.1 Contribución del ala 


La contribución del ala de la aeronave a la estabilidad direccional depende de dos factores, 
estos son, el ángulo de diedro I' y el ángulo de flecha A. Esto mismo se expresa matemáticamente 
como lo dado por la Expresión Se recuerda que si el ala no tiene diedro y no tiene mucho 
ángulo de flecha, la contribución del ala no es tan significativa p. 262). 


(Cs), E (Cn)... + (Cn), y 62) 
de la Expresión [32] la contribución por el ángulo de diedro T' puede estimarse empleado la Ex- 
presión [33] y la contribución debido a la flecha A puede obtenerse mediante la Expresión [34](U.S.| 
Department of Defense||1977). 


(Cs), ,, = 0.075 CL = -0.001501 [1/rad] (33) 
(Cu). dl A A? PELO 
C2 ——4mA  nrA(A+4cos(A./4)) 1472 8cos(A.ys) A 


(Ca), ,, =0,000985 [1/rad] 
finalmente, la contribución del ala resulta: 


(C,g),, =-0,000516 [1/rad] 


La contribución del diedro será desestabilizante, mientras que la contribución debida a la flecha 
es estabilizante. 
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3.3.2.2 Contribución del fuselaje 


La contribución del fuselaje es desestabilizante, se puede aproximar mediante la Expresión [35] 


S l 
(Ca), = KK | 22) (4) [1/deg] (35) 


en donde el factor Ky es un parámetro de carácter empírico por la interferencia existente entre el 
ala y el fuselaje y puede obtenerse mediante el gráfico de la Figura[37/en el cual se entra teniendo 
en cuenta algunos parámetros geométricos del fuselaje. El factor Kg, es un parámetro empírico 
que es función del número de Reynolds basado en la longitud total del fuselaje /¿ y fue obtenido 
mediante la Figura|38| Finalmente la contribución del fuselaje resulta: 


(Cs), = 0, 161436 [1/rad] 


Sa,s: Body Side Area 


Dimar: Maximum Body Width 
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Figura 37: Obtención del factor Ky 
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Ri, fusclage x 1076 


Figura 38: Obtención del factor Kp; 


3.3.2.3 Contribución del estabilizador vertical 


El empenaje vertical, de forma análoga a lo que sucede en la estabilidad estática longitudi- 


nal, aporta la principal contribución para la estabilidad estática direccional de la aeronave. Dicha 
contribución puede ser estimada aplicando la Expresión 


Car OC, dr 
78 = Vit agr = Var + 7) = 0,460113 [1/rad] (36) 


en donde DATCOM se introduce una formula empírica que permite obtener la influencia combi- 
nada del ala-fuselaje en el sidewash según: 


(Crg)r = 


dr e ca 
1 +=] = 0,724 + 3,06————- + 0,4 0,009A = 1,306017 
0 3) y 1 +cosA< nas il 
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3.3.2.4 Contribución de la planta propulsora 


El momento de guiñada producido por la fuerza normal en la hélice es calculado análogamente 
que el incremento de momento de cabeceo en la planta propulsora visto para la estabilidad estática 
longitudinal, en esta situación se tiene: 


Caja ad 0,009174 [1/rad] 
a E TE RS is 
e 
en donde el termino -- es el obtenido en la Figura|25] 


3.3.2.5 Estabilidad Estática Direccional Total - Efecto Veleta 


La derivada C;,g total será la suma de cada una de las contribuciones desarrolladas y calculadas 
anteriormente, esto resulta lo dado en la Expresión [37] 


>| Cs =0,288988 [1/rad] (37) 


3.3.3. Potencia del timón 


La principal superficie de control direccional es el rudder y su efectividad se cuantifica me- 
diante el parámetro C,,5,, el cual se conoce como la potencia del rudder o del timón. Este parámetro 
es igual al coeficiente de momento de guiñada por unidad de deflexión del rudder. Se puede deter- 
minar mediante la Expresión [38] 


>| Cuór ==9 Va, = 0,2184 [1/rad] (38) 


3.3.4. Viento cruzado admisible en el despegue 


Durante el despegue, la presencia de viento cruzado puede generar un deslizamiento lateral en 
el avión, ya que el vector de velocidad resultante se aparta del plano de simetría de la aeronave. En 
casos donde la estabilidad direccional es positiva (C,g > 0), la aeronave tiende a alinearse con la 
dirección del viento para corregir el deslizamiento. Sin embargo, esto puede plantear problemas de 
seguridad durante la carrera de despegue. Para contrarrestar este efecto, se requiere que el timón 
genere un momento de guiñada, esto se esquematiza en la Figura 

Buscamos estimar la velocidad de viento cruzado máxima Vo max posible, teniendo en cuenta 
la deflexión máxima del rudder 0, max y €l C,g a mando fijo. La deflexión del rudder para generar 
el ángulo de deslizamiento máximo posible se encuentra dada en la Expresión [39] (Ver 


2004; p. 281)). 


Cap Braz DOnmár Lañs 
Ona: = 5 Lía > Prez = Pe o = 15, 87 (39) 


Cas Ig Cr 
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Obtenido Bmax, es posible calcular la velocidad del viento cruzado máximo mediante la Ex- 
presión Teniendo en cuenta el valor de la velocidad de despegue calculado en la sección de 
Performance (Ver Sección P.6.3), el valor del viento cruzado máximo resulta lo indicado en la 
Expresión 


Vnmáz 
Dias = tan”! - ) > ado = tan nar) Viof (40) 
Vr 
Pl Venmax= 14,47 [m/s] = 28,12 [nudos] (41) 


PISTA 


E] 
Il 


Figura 39: Diagrama esquemático del viento cruzado durante la maniobra de despegue. Con defle- 
xión del timón nula y con deflexión del timón un ángulo ó,.. 
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3.3.5. Vuelo con un motor inoperativo 


En aviones multimotor, la falla parcial de uno da lugar a una situación de potencia asimétrica 
que puede generar un momento de guiñada significativo, como se muestra esquemáticamente en la 
Figura|40] Si este momento de guiñada no es contrarrestado por el timón, la aeronave desarrollará 
un deslizamiento lateral, el timón debe estar diseñado para tener suficiente autoridad de control 
para mantener el vuelo controlado en tales condiciones. 


| 
y E 
I YE E a 
ll Vb j 
= 


Figura 40: Efectos de potencia asimétrica 


El propósito de esta sección es verificar que el timón pueda contrarrestar esta condición de 
potencia asimétrica en dos condiciones de vuelo, condición de crucero y condición de aterrizaje. 
3.3.5.1 Condición de crucero 

El momento de guiñada producido por el empuje asimétrico se encuentra dado por la siguiente: 

Nr =T-d 


en donde el valor del empuje 7 es obtenido a partir de la curva polar en condición de crucero (Ver 
Sección [6. 1) y d es la distancia del eje del motor al plano de simetría de la aeronave. El valor del 
momento de guiñada resulta: 


Nr = 3306,15 - 3,10 = 10249, 08 [Nm] 
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El momento de guiñada producido por la deflexión del rudder es: 


1 
N, = zPVS DCnórór 
Por condición de equilibrio: 
N7+N,=0 
Finalmente, la deflexión del rudder requerida para mantener el equilibrio se puede encontrar 


mediante la Expresión [42] 


—2N7T 


Li 42 
PS OCr ) 


r 


3.3.5.2 Condición de aterrizaje 
Para la condición de aterrizaje el procedimiento es análogo a la condición de crucero, se utiliza 
la curva polar de aterrizaje para la determinación del empuje. El momento de guiñada resulta: 
N7 = 4741 - 3,10 = 14696, 85 [Nm] 


Finalmente teniendo en cuenta la condición de equilibrio, se determina el ángulo de deflexión 
del rudder requerido para mantener el equilibrio. La deflexión requerida resulta: 


 —2N 
y PaÑaS Cno 


= 4,059 


r 


3.3.6. Discusión 


En esta sección se determinó la estabilidad estática direccional C,,, de la aeronave a mando fijo 
siguiendo la metodología propuesta en la bibliografía 
[1977), se pudo verificar la condición de estabilidad direccional la cual fue descripta al inicio del 
tema (C,, > 0). Se puede concluir, teniendo en cuenta los resultados obtenidos que la aeronave es 
direccionalmente estable. 

Se verificaron las condiciones de viento cruzado admisible en el despegue, los resultados obte- 
nidos arrojan que la aeronave puede tolerar hasta 14.47 [m/s] (equivalente a 28.12 [kn]) de viento 
cruzado; el fabricante de la aeronave informa que el valor máximo admisible es de 15.43 [m/s] 
(equivalente a 30 [kn]). Se considera que el cálculo realizado arroja resultados satisfactorios. 

Finalmente otro requerimiento verificado es el vuelo con un motor inoperativo, se determinó la 
deflexión del rudder requerida para poder mantener a la aeronave trimada direccionalmente. Los 
ángulos de deflexión requeridos obtenidos están dentro de las deflexiones posibles para el timón, 
resulta algo llamativo que los valores de dichos ángulos sean relativamente pequeños. 
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3.4. Estabilidad Estática Lateral 


3.4.1. Introducción 


En esta sección, se lleva a cabo el análisis de la estabilidad lateral de la aeronave asignada al 
grupo desde una perspectiva estática. El principal control sobre el momento de rolido lo propor- 
ciona la familia de alerones, ubicados cerca de las puntas de las alas. Al deflecionarse, aumentan 
la curvatura del ala, afectando la distribución de presiones y generando una variación en la susten- 
tación que resulta en un momento de rolido. En la Figura 41] se presenta la distribución de carga 


en relación con la velocidad de rolido (Perkins y Hagej|1949). 


L 


D (radians per a 


A 
AS A 


Figura 41: Distribución de carga debido a la velocidad de rolido 


3.4.2. Criterio de estabilidad 


Para que la aeronave sea estáticamente estable alrededor del eje x debe desarrollar un mo- 
mento de rolido que restablezca el ángulo de inclinación original. El momento de rolido que se 
produce debido al deslizamiento se conoce como efecto diedro y el criterio para un efecto diedro 
estabilizante es: 


Cr <0 


3.4.2.1 Contribución Ala - Fuselaje 


La contribución del ala de la aeronave para la estabilidad lateral depende principalmente del 
ángulo de diedro I” y el ángulo de flecha A, mientras que la contribución del fuselaje se ve reflejado 
por diferentes factores como se detalla en la Expresión [43|del método DATCOM: 
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E AC 
er Ko + == 


+ (AC;,) (43) 


Zw 


Donde: 
C£g 


La contribución del ala debida a la flecha (5) se determina mediante la gráfica de la Figura 
c/2 


SUBSONIC SPEEDS 


A¿¡2 (deg) 
20 0 20 40 60 80 


(3 
CL Mera * 


(per deg) 


Aga (deg) 


Figura 42: Contribución del ala debida a la flecha 


(E) =0 
Cl Aoe/2 


E 
La contribución debida al alargamiento (E) se determina empleando la gráfica de la Figura 
A 


C 
(E) = -0,001 
C A 


L 


Ahora, la contribución del ala debida al diedro a se obtiene de la gráfica de la Figura[44] 


Ct; 
—= == 2 
T 0,00025 


El factor K yr se determina mediante la gráfica de la Figura[45] 


Kur = 1,04 


Se obtienen las contribuciones finales de la interferencia ala-fuselaje: 


AC £, 


d 2 
= 0,0005 VA (5) = -0,00003 (44) 
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Aspect Ratio A 


Figura 43: Contribución del ala debida al alargamiento 


-.0003 y-— 


Sa -.0002 


(per deg?) 


Figura 44: Contribución del ala debida al diedro 


SUBSONIC SPEEDS 


Figura 45: Correcciones por compresibilidad debida a diedro 
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Donde d es el ancho máximo del fuselaje. 


1,2 VA (2,1 (2d 
co), E) oso ás 


Donde z,, es la posición del ala relativa al eje x. 


Reemplazando todas las contribuciones y correcciones obtenidas en la Expresión|43] se obtiene 
la contribución total del ala-fuselaje al efecto diedro, quedando: 


(c o = -0,1488 [1 /rad] 


3.4.2.2 Contribución del empenaje 
La contribución del estabilizador vertical puede ser expresada como la Ecuación [46] 


(Ca) A 17N (46) 


Donde la distancia z, se indica en la gráfica de la Figura 


Figura 46: Posición del estabilizador vertical 


Reemplazando los datos correspondientes del estabilizador vertical en la Expresión[46] se ob- 
tiene la contribución del estabilizador vertical: 


(C:,), = -0,1022 [1/rad] (47) 


3.4.2.3 Estabilidad Estática Lateral Total - Efecto Diedro 


La derivada Cy, total será la suma de las contribuciones calculadas en las secciones anteriores, 
esto resulta lo dado en la Expresión|48] 


>| C¿,=-0,2510 [1/rad] (48) 
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3.4.3. Potencia de alerones 


Para obtener la Potencia de alerones se hace uso de la Expresión [49] 


Ct, = f cydy (49) 
k 


Donde k, y k, son las posiciones a lo largo de la envergadura de las estaciones límites del 
alerón,como se aprecian en la gráfica de la Figura[47|y 7, €s el factor de efectividad de los alerones. 


p b/2 


k1b/2 K2b/2 


Figura 47: Parámetros del alerón 


El resultado se obtiene mediante la integración numérica aplicando la Regla de Simpson ha- 
ciendo uso de Python, reemplazando con los datos correspondiente del ala, se obtiene: 


»| Cr, =0,3272 [1/rad] (50) 
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3.4.4. Velocidad de rolido estacionaria 


La velocidad de rolido estacionaria se alcanza cuando el momento debido a la deflexión de 
alerones iguala al momento de amortiguamiento debido a la propia velocidad de rolido. 
3.4.4.1 Método NACA 

Haciendo uso del método NACA propuesto en el Report NACA TR 635, se tiene la Ecuación|51] 


para obtener la velocidad de rolido estacionaria : 
pb Cr, TÓLK 
2V  T1146C,, 


(5D) 


Donde los parámetros E Cr, 7 y K, se obtienen de las Figuras [50] y [51|respectiva- 
mente. 


2 A .6 .8 1,0 

Extent of dihedral, 9% 9/2 Ñ 

TFicurE 9-3. Chart for estimating C1y due to wing dihedral (TP), From TR 635, 

“Theoretical Stability and Control Characteristics of Wings with Various Amounts 
of Taper and Twist,” hy Pearson and Jones, 


: a E 
Figura 48: Parámetro — 
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Aspect ratio 


Figura 49: Parámetro C£, 


Ca/Cw 


Figura 50: Efectividad del alerón 


Figura 51: Corrección por deflexión del alerón 
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Como se aprecia en las imágenes, se obtienen los parámetros correspondientes: 


C 
o 
T 


C, = 0,62 
T=0,55 
k=0,76 


Tomando en cuenta que la deflexión total del alerón ó, es de 41” y reemplazando todos los 
parámetros obtenidos en la Expresión|[51) se llega a la velocidad de rolido estacionaria: 


pe = 0,06512 (52) 


2V 


Despejando de la ecuación [52] la velocidad de de rolido dimensional resulta: 
p=0,5713 [rad/s] = 32,73 [deg / s] 


3.4.4.2 Método de integración por franjas 


También se puede obtener la velocidad de rolido estacionaria haciendo uso de la Expresión|[53] 


k23 
pb Tab dz cy dy 


= (53) 
b 
2V Z le cy? dy 
La cual puede ser resuelta mediante la teoría de integración por franjas a partir de: 
pb ÓTa (BTRUS — 4) + (1 — TRIOS — ki) (54) 
Y SS 3+TR 
Donde 
C, 
TR = = 
Cr 


Y k¡ y k, son las posiciones a lo largo de la envergadura de las estaciones límites del alerón, 
como se puede ver en la Figura[47] 
Se llega entonces a la velocidad de rolido estacionaria por el método de integración por franjas: 


(55) 
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3.4.5. Discusión 


En este apartado se logró determinar la estabilidad estática lateral de la aeronave C 1, siguiendo 
la metodología DATCOM. Se comprobó el cumplimiento de la condición de estabilidad lateral 
presentada C ,, < 0, concluyendo así que la aeronave es lateralmente estable. 

Para la determinación de la estabilidad estática lateral C.,, se obtubieron las contribuciones 
tanto del ala como del empenaje vertical. 

Así mismo se obtuvo la potencia de alerones C,,, y la Velocidad de rolido estacionaria adi- 
mensional. Esta última fue determinada mediante dos métodos, el Método NACA y el Método de 
integración por franjas. Se observó una discrepancia en los resultados obtenidos entre ambos mé- 
todos. Es importante destacar que el Método de Integración por Franjas, al no considerar efectos 
tridimensionales, tiende a proporcionar valores mayores, lo cual es normal. Por lo que se considera 
que los resultados son acertados. 
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4. Estabilidad Dinámica 


4.1. Introducción 


En las siguientes secciones, se abordará el análisis de la estabilidad dinámica de la aeronave 
asignada al grupo, el Short 360. Para llevar a cabo este análisis, se empleó la metodología propuesta 


por Etkin (Etkin y Reid|[1996) y para la estimación de las derivadas de estabilidad adimensionales 
se emplearon las metodologías propuestas por dicho autor y las metodologías presentada en el 


DATCOM (U.S. Department of Defense| 1977). 


4.2. Estabilidad Dinámica Longitudinal 
4.2.1. Introducción 


En esta sección, se aborda el análisis de la respuesta dinámica longitudinal del Short 360, 
la aeronave asignada al grupo, para pequeñas perturbaciones. Siguiendo las metodologías reco- 
mendadas por la cátedra y las referencias bibliográficas 
[Defensej [1977), se estimaron las derivadas adimensionales longitudinales. Posteriormente, se cal- 
cularon las derivadas dimensionales y se construyó la matriz A, que permitió resolver el sistema 
dinámico. Los valores y vectores propios asociados a los dos modos longitudinales, el fugoide y 
el corto período, fueron finalmente obtenidos. Los datos empleados para resolver esta sección se 
encuentran dados en la Tabla 


Tabla 11: Datos necesarios - Estabilidad Dinámica Longitudinal 


Parámetro Símbolo Valor Unidad 
Coeficiente de sustentación en crucero Cr 0.4866 Adim. 
Coeficiente de resistencia en crucero Cp 0.0315 Adim. 
Pendiente de la configuración completa a 6.8687 1/rad 
Pendiente de sustetentación del E.H. a; 4.3523 1/rad 
Posición del C.G. en crucero h 0.2650 Adim. 
Densidad en el crucero Yo) 0.9044 — kg/m' 
Punto neutro de mando fijo con potencia h 0.5635 Adim. 
Velocidad de crucero U0 100.0000 m/s 
Mach en crucero M 0.3045 Adim. 
Volumen de cola Vi 1.2182 Adim. 
Distancia entre C.A. empenaje horizontal y C.G. del avión l, 10.5377 m 
Cuerda Aerodinámica Media Cc 1.9900 m 
Gradiente de Downwash de/da 0.3640 Adim. 
Superficie alar y 42.7300 m? 
Momento de inercia en y L, 121553 kgm? 
Ángulo de actitud en crucero 00 0.0000 > 
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4.2.2. Derivadas de estabilidad adimensionales 


Se estimaron siguiendo la metodología desarrollada en la bibliografía de referencia (Ver 
1996) p. 155)). Las fórmulas empleadas se encuentran sintetizadas en la Tabla 


Tabla 12: Derivadas longitudinales adimensionales 


Ex C, Cn 
0 ME) AU CI) ME) PUC 
Q C.- Cho, =(C;, + Cp) -a(h, — h) 
d 0 24 VH E 24 Vni% 
q 0 —2a,Vy -2a,Vy 


Considerando como despreciables las derivadas C.,,, Cx, y Cm,» y asumiendo que los efectos 
aeroelásticos son despreciables, se obtuvieron los resultados presentados en la Tabla En este 


análisis, para una aeronave turboprop, Cy, se aproxima como —3C7,, donde Cr, = Cp +C;,, sin(0p)) 


- Y 
y Emp = Pas * 


Tabla 13: Resultados de las derivadas adimensionales longitudinales 


E C. Cn 
u -0,0946 -0,0497 0,0000 
a 0,2480  -6,9002  -2,0505 
0,0000  -3,8773  -20,5316 
q 0,0000 -10,6519 -56,4055 
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4.2.3. Derivadas de estabilidad dimensionales 


Posteriormente, al considerar la condición de vuelo crucero y los resultados de las derivadas 
adimensionales de estabilidad longitudinal, se determinaron las derivadas de estabilidad dimensio- 
nales siguiendo la formulación establecida en la Tabla [14] 


Tabla 14: Derivadas Dimensionales Longitudinales 


X Z M 

u — puyS Cy, sin Oy + 3 ppuoS Cy,  —puoS Cy, COS Op + 3puoS Co,  ¿puyES Cm, 
w ¿puoS Cae ¿puoS Cas ¿PUYTS Ca 
q ¿pupas Cs, ¿puyasS Co, ¿puy??S Cm 
w ¿pas Ca, ¿pas Co, ¿pS Con 


Los resultados obtenidos se presentan de forma similar a la sección anterior en la Tabla|15] 


Tabla 15: Resultados de las Derivadas Dimensionales Longitudinales 


X [NM] Z[N] M [Nm] 

u [m/s] —182,7909 —1976,4991 0,0000 
w[m/s]  479,1981  -—13332,9156  -—7884,5389 
q [rad/s] 0,0000 -20479,2295 -215804,7807 
w [m/s?] 0,0000 -74,5446 -785,5293 


4.2.4. Modos longitudinales 


El sistema de ecuaciones de pequeñas perturbaciones para los modos longitudinales es el que 
se encuentra dado por la Expresión [56] 


Xu Xy - 
ú m m 0 g cos (00)| y y, 
. Zu Zo Z¿+muo gm sin (00) 
w 7 —Z yy +m —Z yy + —Z y +m Zim w ( 5 6) 
Gl | MuZo-m)-MiZ, — Mu(Zim)-MyZ Ma(Zy 1) -My(Z¿+muo)  Mygmsin(00) q 
1) Ly(Z;—m) Ly(Z—m) Ly(Z im) Ly(Zy—m) 0 
0 0 1 0 
De forma sintética puede ser expresado como se muestra en la Expresión [57] 
x= Ax (57) 


en donde x es el vector de estado y A es la matriz del sistema la cual posee valores constantes. Las 
soluciones de este sistema de ecuaciones diferenciales de primer orden tienen la forma ya conocida 
como: 


x(t) = Xpe* 
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En la solución Xy es el vector propio y 1 es el valor propio del sistema. En síntesis, solo habrán 
soluciones cuando se cumpla lo siguiente: 


det(A—211)=0 


El determinante de la expresión estipulada anteriormente dará como resultado la ecuación 
característica del sistema, la misma será un polinomio en términos de 1 de cuarto orden y sus 
raíces serán los valores propios del sistema. 

Aplicado al caso de estudio que se está tratando actualmente, la matriz A del problema, resulta 
la dada en la Expresión|[58| 


-0,019066 0,049983  0,000000  -—9,806650 

—0,204568 -—1,379960 97,108858  0,000000 
0,001322  —0,055947 -—2,402957 0,000000 
0,000000  0,000000  1,000000  0,000000 


Su ecuación característica resulta: 


A= (58) 


A + 3,8020? + 8,8313 + 0,19792 + 0,1301 = 0,0000 


Calculando las raíces de la ecuación característica, se obtuvieron los siguientes valores propios: 


Modo Fugoide 112 = -0,0081 + 0,1216; 
Modo Corto Período 234 = —1,8929 + 2,2742 ¡ 


Vemos que los modos naturales son dos oscilaciones amortiguadas, una de período largo y 
ligeramente amortiguada (el Modo Fugoide) y la otra de período corto y fuertemente amortiguada 
(Modo de Corto Período). Los vectores propios asociados a los valores propios obtenidos son los 
siguientes: 


0,999058 + 0,01 0,013980 + 0,002377¡ 

X= 0,041535 + 0,0006701 Xi 0,999579 + 0,01 
127 |-0,001519 + 6,149161£ — 05i 34 — |-0,005250 + 0,023404i 
0,001330 + 0,0123961 0,007215 + 0,0036961 


Obtenidos los vectores propios anteriores, los normalizamos de forma tal que O = 1. Para ello 
dividimos cada elemento de cada vector propio por el valor de 6. El resultado de esto último lo 
presentamos a continuación: 


8,5481 + 79,6778i 1,6686 + 0,5254i 
ia 0,4088 + 3,3068¡ opio 109,7428 + 56,2245¡ 
12 | -0,0081 + 0,1216 34 — | 1,8929 + 2,2742 


1,0000 + Oi 1,0000 + Oi 


Finalmente se expresan los vectores propios normalizados en su forma adimensional, para 
realizar dicha adimensionalización se sigue lo siguiente: 


Desimone - Álvarez Granda 85 de 


Trabajo Integrador Grupal - Short 360 


" u 
Uu=— 
U0 
== 
uN 
G= q 
2 -upy/c 


Los vectores propios normalizados adimensionales finalmente, son: 


0,0855 + 0,7968 0,0167 + 0,0053i 

ye, | 0,0041 + 0,0331 xa. — | 1,0974 + 0,5622i 

12 [0,0001 + 0,0012í 34 — 120,0188 + 0,0226i 
1,0000 + Oi 1,0000 + Oi 


Los vectores y valores propios son indeterminados en magnitud, por lo tanto, podemos mul- 
tiplicarlos por cualquier número (real o complejo). En este caso, elegimos que la componente 
asociada a O valga 1. De esta manera, las perturbaciones de velocidad u, ángulo de ataque w, y 
velocidad angular de cabeceo q mostrarán el ángulo de fase respecto al ángulo de actitud del avión 
y la magnitud relativa a un cambio unitario en el ángulo de actitud. Esto resulta útil porque pue- 
de visualizarse en un diagrama de fases, como se muestra en la Figura [52] los vectores propios 
pasados a forma polar se presentan en la Tabla[16] 


Tabla 16: Vectores propios en forma polar 


Modo Fugoide Modo de Corto Período 
Módulo Fase Módulo Fase 
ú 0.8014 -96.1? 0.0175 17.59 
w 0.0333 82.9” 1,2331 2 
gq 0.0012  -93,8* 0.0294 129.8" 
68 1.0000 0.02 1.0000 0.09 


Los diagramas fasoriales de la Figura[52]permiten corroborar la teoría. En el modo fugoide, se 
considera un movimiento en el que la velocidad de cabeceo q y el ángulo de ataque « = Y» son muy 
pequeños (y de no ser por la amplificación no serían percibidos en el diagrama), pero ú y O están 
presentes con una magnitud significativa. En cuanto al modo de corto período, las variaciones en 
velocidad son despreciables, mientras que el ángulo de ataque a = Y y 0 están presentes con una 
magnitud significativa. 

Finalmente empleando la expresión general de soluciones se construyeron los gráficos para las 
distintas variables de interés longitudinales, esto se presenta en las Figuras [53a] y [53b] 
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Diagrama de fases modo fugoide 


9 - Amplificado (No es visible) 
- Amplificado (No es visible) 


Diagrama de fases del modo corto período 


ú - Amplificado (No es visible) 


- Amplificado (No es visible) 


Figura 52: Diagramas fasoriales de los modos fugoide y corto período 
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O 100 200 300 400 500 600 O 100 200 300 400 500 600 
t — [s] t — [s] 


O 100 200 300 400 500 600 O 100 200 300 400 500 600 
t — [s] t — [s] 


00. 05 10 15 20 25 30 35 00. 05 10 15 20 25 30 35 
t — [s] t — [s] 


(b) Modo de Corto Período 


Figura 53: Modos longitudinales 
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Finalmente se determinaron los períodos, los amortiguamientos y tiempos de reducción a la 
mitad, los mismos se presentan en la Tabla[17] 


Tabla 17: Períodos, amortiguamientos , tiempos de reducción a la mitad y número de ciclos para 
reducción a la mitad. 


Modo Período Tiempo de reducción a la mitad Amortiguamiento 
[s] [s] [-] 
Fugoide 51.65 85.68 0.0664 
Corto Período 2.76 0.37 0.6397 


4.2.5. Modos longitudinales de orden reducido 


= Modelo de Lanchester: En el modelo de Lanchester se desprecia el cambio de ángulo de 
ataque, y se cosidera que las fuerzas aerodinámicas no realizan trabajo en la dirección de la 
trayectoria. Esto se traduce en la conservación de la energía mecánica 


E =1/2mV? — mgz, 


De donde se obtiene la velocidad en términos de la velocidad de vuelo y su cambio de 
velocidad debido a la variaci?on de energía potencial gravitatoria. 


Dz 
V" =u, +28Z. 
Relacionando L y W ,y reemplazando en las ecuaciones de movimiento se obtiene: 


W - Lcos(y) = W - L = mz, 


Donde vemos que la diferencia entre W y L es debida al cambio de velocidad por variación 
de la energía potencial gravitatoria y es proporcional a z,, por lo que el sistema pasa a ser 
análogo a un sistema masa resorte sin amortiguamiento con rigidez de resorte equivalente 
igual ak = C,.pgS y queda expresado como: 


mz. +kzg =0 


m 
AN 
"Nk 


= Modelo de orden reducido para modo fugoide: En este modelo el momento de cabeceo 
se asume muy pequeño y que el avión se encuentra en una sucesión de equilibrios cuasies- 
táticos. De esta manera tenemos que q = 0 y el sistema se reduce en un grado de libertad, 


y en donde el período T resulta: 
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además teniendo en cuenta que 0, = O se resuelve un problema de autovalores obteniéndose 
una cuadrática: 


AR+BA+C=0 


A = —u0M., 
B= 32M, + (X,M, — M,X,,) 
m 


C=É(Z,M,-M,Z,) 
m 


Además si consideramos que Cm, = 0, se simplifica aún más: 


= Modelo de orden reducido de corto periodo: Se asume que la velocidad se mantiene 
constante y se toma en consideración que m es mucho más grande que Z,, y que mu, es 
mucho mayor que Z, para simplificar las ecuaciones. Se llega a un polinomio característico 
de segundo grado. 


P+BA+C=0 


A | 
B=- z + —UM, + Mo) 
ml, 


1 
o lun, 
L m 


y 


Resolviendo los modelos de orden reducido se obtuvieron los resultados presentados en la 
Tablas[18] y [19] Los valores propios obtenidos resultaron: 


Modo Fugoide 212 = -0,0095 + 0,1419 
Modo Corto Período 234 = —1,9062 + 2,3070 ¡ 
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Tabla 18: Modelo de orden reducido para Modo Fugoide 


Parámetro Ecuaciones generales Lanchester  q=0 
Frecuencia amortiguada [rad/s] 0,1216 0,1387 0,1419 
Período [s] 51,6519 45,3038  44,2891 
Relación de amortiguamiento 0,0664 - 0,0670 
Tiempo 1/2 [s] 85,6806 - 72,7103 


Tabla 19: Modelo de orden reducido para Modo de Corto Período 


Parámetro Ecuaciones generales Modelo de orden reducido 
Frecuencia amortiguada [rad/s] 2,2742 2,3070 
Período [s] 2,7628 2,7236 
Relación de amortiguamiento 0,6397 0,6370 
Tiempo 1/2 [s] 0,3662 0,3636 


4.2.6. Discusión 


Analizando los resultados obtenidos se puede observar que los modos obtenidos se ajustan a lo 
esperado; uno de ellos presenta un período relativamente extenso y un amortiguamiento reducido, 
correspondiente al Modo Fugoide, mientras que el otro muestra un mayor amortiguamiento y un 
período corto, característico del Modo de Corto Periodo. Demostrando así que ante una perturba- 
ción este último se estabilizará mucho más rápido que el fugoide. 

En ambos casos, se evidencia estabilidad dinámica, ya que la parte real de los valores propios 
asociados es negativa y como se pudo apreciar en las gráficas presentadas existe una evidente 
convergencia. 

Comparando los resultados obtenidos con los modelos de orden reducido (Tablas [18] y [19] se 
percibe una destacada similitud numérica, enfatizando la significativa proximidad entre los valo- 
res obtenidos. Esta observación respalda y confirma la consistencia de los resultados calculados, 
subrayando así la confiabilidad de los mismos. 

Las formulaciones matemáticas más sencillas de los modos de orden reducido nos permitieron 
obtener valiosas conclusiones. Al observar en particular el modo de corto período, se evidencia su 
fuerte influencia por M, (C,,). A medida que este valor aumenta, la parte real, que se corresponde 
con el amortiguamiento, adquiere mayor importancia. La derivada c,,, está vinculada al volumen 
de la cola Vy, y dado que el volumen de cola de la aeronave es significativo (superior a 1.20 P, 
es lógico que se manifieste un valor elevado de M,. Además, la parte imaginaria (frecuencia) es 
considerable, posiblemente debido a que M,, (relacionado con c,, ) también tiene un valor alto. En 
este caso, a medida que la derivada c,,, aumenta, la parte imaginaria adquiere mayor magnitud. Es 
importante destacar que nuestro alto margen estático contribuye a obtener un c,,, elevado. 


*Dicho valor de volumen de cola es característico de aeronaves tipo regionales (Turboprop Regional), según indica 
el autor Schaufele en su libro, valores típicos para aeronaves de este tipo están comprendidos entre 0.83 y 1.47 


p- 147, 
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4.3. Estabilidad Dinámica transversal 


4.3.1. Introducción 


En la sección presente se aborda el análisis de la respuesta dinámica transversal para pequeñas 
perturbaciones de la aeronave Short 360. Siguiendo las metodologías recomendadas por la cátedra 
y las referencias bibliográficas (Etkin y Reid|[1996;|Roskamj[1973), se determinaron las derivadas 
adimensionales laterales direccionales. Seguidamente se calcularon las derivadas dimensionales y 
se construyó la matriz A, que permitió resolver el sistema dinámico. Finalmente se obtuvieron los 
valores y vectores propios asociados a los modos laterales direccionales, el modo de rolido puro, 
el modo de espiral (Spiral), el modo de balanceo del holandés (Dutch roll). Los datos necesarios 
para la resolución de esta sección se encuentran en la Tabla|[20] 


Tabla 20: Datos necesarios - Estabilidad Dinámica Transversal 


Parámetro Símbolo Valor Unidad 
Ángulo de diedro P 3.0000 Á 
Coeficiente de sustentación en crucero Cr 0.4866 Adim. 
Coeficiente de resistencia en crucero Cb 0.0315 Adim. 
Densidad en el crucero p 0.9044  kg/m* 
Velocidad de crucero U0 100.0000 m/s 
Mach en crucero M 0.3000 Adim. 
Pendiente de sustentación E. V. a, 4.2288 1/rad 
Superficie del E. V. Sr 7.7300 m? 
Pendiente de la configuración completa a 6.8687 1/rad 
Alargamiento A, 12.1657  Adim. 
Distancia l, Le 10.7760 m 
Distancia z,, > 2.8140 m 
Ángulo de incidencia de crucero Q 4.0606 Ñ 
Momento de inercia en x Le 104785 kgm? 
Momento de inercia en y L, 121553 kgm? 
Momento de inercia en z L 210172 kgm? 
Envergadura del EH. b, 7.4800 m 
Superficie del E.H. S» 9.8300 m? 


4.3.2. Derivadas adimensionales de estabilidad transversal 


Las derivadas de estabilidad lateral direccional se obtuvieron siguiendo la metodología desa- 


rrollada en la bibliografía de referencia (Ver 1973)). Las fórmulas empleadas se detallan 
a continuación. 


Desimone - Álvarez Granda 92 de [126] 


Trabajo Integrador Grupal - Short 360 


4.3.2.1 Derivadas con respecto a 5 


Cy 


B 
La derivada Cy, total es la suma de las contribuciones correspondientes: 
Cy, = Cry, + Crg, + Cro, (59) 
La contribución del ala se determina a partir de la Expresión[60] 


Cy 


B 


y = 70,0001/1157,3 [1/raa] (60) 
La contribución del fuselaje: 
S 
Cray =2K:(2) [1/rad] (61) 


Donde $ es la sección transversal del fuselaje en xy donde se separa el flujo. xy se obtiene en 
función de x, como se aprecia en la Expresión [62], el punto donde la pendiente E del fuselaje 
tiene mayor valor negativo. 


== 0,378+0,5272 (62) 
l, l, 
Y finalmente, la contribución del empenaje vertical: 
00 S, 
Cy, = Como se 7) al [1/rad] (63) 


Ct; y Cup 


La obtención de las derivadas Cr, y C,, se detallan en los apartados de estabilidad estática 
lateral y direccional correspondientemente. 


4.3.2.2 Derivadas con respecto a la velocidad de rolido p 


Cy, 
La derivada Cy, se obtiene mediante la Expresión|64] cuya contribución más significativa es la 
del empenaje vertical. 
zycos(a) — l,sinlo 
O ORT 
Donde la posición y distancia z, y l, se determina como lo muestra la Gráfica [54] 
Ct, 


Cs (64) 


v 


La derivada C 7, total, se obtiene de la suma de las contribuciones del ala-fuselaje, empenaje 
horizontal y vertical. 
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Figura 54: Parámetros de posición del empenaje vertical 


BODY Xx-AxiS 


Xa 


AIRPLANE C.6, 


La contribución del ala-fuselaje se determina con la Expresión: 
BCr 
a = - 


K 


B 


K 


Donde: 


pi E 
El parámetro (E =p 


K 


) se obtiene de la Gráfica 


(d) A = 1.0 


Ef] 


, 
¿2 


ll 


: 


A, (deg) = tan» 
y (des) = tan ” (tama, /,/8) 


ea ) 


Figura 55: Parámetro ( 


(65) 


El parámetro x es la relación entre la pendiente de sustentación 2D(del perfil) a la teórica: 


(04 


277 


C) 
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El parámetro f£ de Prandtl Glauert de corrección por compresibilidad, se definde como: 


B= xy/1 - M*cos(Agos) 


La contribución del empenaje horizontal se determina con la siguiente expresión; 


2 
Ce, =0,5(C:,,) ad (2) (66) 


Donde la derivada (C £ru) se obtiene con el mismo procedimiento de la contribución del ala- 
fuselaje que se detalló en el apartado anterior, cambiando los datos del ala con los del empenaje 


horizontal. 
La contribución del empenaje vertical se calcula con la ecuación: 


Cn (Or) (2) (67) 
C,, 


La contribución más significativa para determinar la derivada C 
con la Expresión [68] 


es la del ala, que se calcula 


Mp 


Ca 
Crow a C£,ytan(a) mn [Cant ES (E) Je. (68) 
CL C¿=0 


Donde: 


A A 5 y Copy 
C, corresponde al coeficiente de sustentación en condición crucero y el parámetro | ¿2 se 
C,=0 


obtiene con la expresión siguiente: 


non [ÁS AAA (69) 


Xcg—Xcaw anAo, an? 
(E) ¡ AR +6(AR + costAnas) | o ( "Aozs al As ) 
Ci=0 6 AR + 4cos(Ao25) 


4.3.2.3 Derivadas con respecto a la velocidad de guiñada r 


Cy, 
La principal contribución para obtener la derivada Cy, la realiza el empenaje vertical, con la 
expresión siguiente: 


Cy =Cy =-2 l,cos(a) + z,sinta) 


“Cr, o Cra, (70) 


v 


Ca in 


La derivada C,, total se obtiene de la suma de las contribuciones del empenaje vertical y el ala. 
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La contribución del empenaje vertical se calcula con la siguiente expresión: 


_ ,cos(a) + zysin(a)) 


Ch, = LC, (711) 
Y la contribución del ala, se determina con la ecuación presentada a continución: 
Cr, = [Lu] 0? + [E] e (72) 
Mw de L Cp, Do 


z . Ca Car . y . 
Donde los términos 2] y (Ez) se obtienen de las Gráficas 5ó|y 57 correspondientemente, 
£ 0 


o OT x is distance from the cg. to the a.c., positive for the a.c, aft of the c.g. 
€ is the mean aerodynamic chord. 


RATIO! Y 
A 
e CIEN (GUS ISI: IEEE: TEE MG OO UE DIO DN E O 
C.,/C,? (per rad) 
, A ¡a 
Figura 56: Parámetro | >" 
L 


Ct, 
De la misma manera que en el inciso anterior, la derivada C 7, es la suma de las contribuciones 
del empenaje vertical y el ala. 
La contribución del empenaje vertidal se determjina como: 


(costa) + z,sinltoa)YMz,costa) — l,sinla)) 


Ct E 2 pb? Ya, 


(13) 
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ASPECT RATIO, A 
4 6 


Figura 57: Parámetro 22) 
0 


La contribución del ala se obtiene de la expresión siguiente: 


C ARseníA 
ca) cp tina y 
Ch C,=0 12(AR + Acos(Ap25) 


(14) 


Una vez calculadas todas las contribuciones de cada derivada, los resultados obtenidos se los 
presenta en la Tabla[21] 


Tabla 21: Resultados de las derivadas adimensionales transversales 


o or EN 
fB -1,0358 -0,2510 0,2890 
p -0,4210 -0,5646 -0,0404 
r 0,6719 0,1965  -0,3305 
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4.3.3. Derivadas dimensionales de estabilidad transversal 


Considerando la condición de vuelo crucero y los resultados previamente obtenidos de las 
derivadas adimensionales de estabilidad transversal, se determinan las derivadas dimensionales de 
estabilidad transversal, haciendo uso de las expresiones detalladas en la Tabla [22] 


Tabla 22: Derivadas Dimensionales Transversales 


Y L N 
5puyS Cro ZpuobS Cr 5puobS Cno 
ZpusbS Co Zpuyb?S Ca, ¿puyb?S 20m 


r ¿puobS Es ¿puob?S C,, ¿pupb?S Ca, 


=S = 


Los resultados obtenidos se presentan en la Tabla[23] 


Tabla 23: Resultados de las Derivadas Dimensionales Transversales 


Y [N] L[Nm] N [Nm] 
v [m/s] -2001,4542 11057, 8837 12731,4570 
p[rad/s] -—9273,4233  -—282273,6729  -20315,2268 
r[rad/s]  14800,1624 98673,2413  -—165967,1161 


4.3.4. Modos transversales 


El sistema de ecuaciones de pequeñas perturbaciones para los modos transversales se encuentra 
dado por la Expresión[75] 


v » Lo Ll Uy  gcos(0p)|[v 
7 , L, , Ir , 
p == E il EN, + LN, 7 + EN, 0 pe (75) 
HO IES EL 2+EL, Esti, 0 ||" 
. E rá ó E AS 
10) 0 1 tan (Op) 0 10) 
La misma se puede expresar como: 
x= Ax 
Donde la Matriz A para el caso de estudio actual se muestra en la Expresión [76] 
0,2088 —0,9673 -—98,4563 9,8067 
e 0,1033 -—2,6956 0,9134 0,0000 (16) 
| 0,0588 —0,1439 -—0,7742 0,0000 
0,0000 1,0000 0,0000 0,0000 
La ecuación característica se obtiene de la Expresión [77 
det(A — AI) =0 (77) 
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Resolviendo la ecuación característica se obtienen cuatro autovalores, un par de valores com- 
plejos conjugados y dos reales, que conforman los modos de la aeronave en el análisis dinámico 
transversal: El par de valores complejos conjugados corresponde al modo denominado Balanceo 
del Holandés, mientras que el valor real más grande corresponde al modo rolido puro y el más 
chico es el correspondiente al modo espiral. Quedando: 


= 1 correspondiente al modo espiral (Spiral) 
= 2, correspondiente al modo rolido puro 


= 13, correspondiente al modo de balanceo del holandés (Dutch roll) 


4.3.4.1 Modo Balanceo del Holandés (Dutch Roll) 


Como se aprecia en la Figura es un modo oscilatorio con un gran acoplamiento entre el 
movimiento de guiñada y de rolido. 


4.3.4.2 Modo Rolido Puro 


Este modo tiene un comportamiento no oscilatorio y está desacoplado de los demás movimien- 
tos como se muestra en la Figura [58c]Se encuentra expresado por un valor propio real y de valor 
absoluto relativamente grande. El amortiguamiento es generalmente grande. 


4.3.4.3 Modo Espiral 


El modo espiral es no oscilatorio y se caracteriza por ser un movimiento lento que involucra 
complicados movimientos de rolido, guiñada y deslizamiento (Ver Figura [58b). Tiene un valor 
propio real y de módulo pequeño. Físicamente el avión desarrolla un giro descendente a medida 
que su nariz cae con respecto a la horizontal. 
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(a) Modo Dutch Roll 


. — Steadily increasing roll angle , ES 

A, v 
—T— S Á Restoring rolling moment 
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Yawing moment | Y B Steadily mcreasing yaw 
due to fin lifi Yo . 


!/ 


Port wing Starboard wing 
Reduction in incidence Increase in incidence 


(b) Modo Espiral (c) Modo Rolido Puro 


Figura 58: Modos transversales 


Resolviendo la Expresión|77) se obtiene: 


244 3,67851* + 8,63241? + 18,52171 + 0,2578 =0 (78) 


Los valores propios que se obtienen a partir de la misma se expresan como: 


Modo Espiral 1, = -0,014011 
Modo Rolido Puro 2, = -2,895677 
Modo Balanceo del Holandés 234 = -0,384404 + 2,4913871 


Los vectores propios asociados a los valores propios obtenidos son los siguientes: 
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0,7729 + 0,01 0,8874 + 0,01 0,9992 + 0,01 
_ |-0,0088 + 0,01 y = [70:4357 + 0,01 x= | 0,0250 + 0,0172i 
11 0,0615 + 0,01 2 —|-0,0050 + 0,01 34 — [0,0029 + 0,0246i 

0,6315 + 0,01 0,1505 + 0,01 (0,0083 + 0,0088 


Una vez obtenidos los valores propios, se procede a normalizarlos respecto a / y a adimencio- 
nalizarlos siguiendo: 


P=v/4w  p=p/Quo/b) += p/Quo/b) 


Los mismos resultan: 


0,0122 + 0,01 0,0590 + 0,01 0,5697 + 0,6036i 
y, = [7 0.0016 + 0,0 y = [70:3301 + 0,01 ye. , — [0,0438 + 0,2840 
[0,0111 +0,0í 2 — |-0,0038 + 0,01 34 — [0,1508 + 0,1794i 
1,0000 + 0,01 1,0000 + 0,01 1,0000 + 0,0000í 


Tomando los vectores propios normalizados adimencionalizados y utilizando la expresión ge- 
neral de solución se graficaron las variables de interés, esto se muestra en las Figuras y 


0.012 0.0000 Lorna 
0.0002 4 
0.010, 
0.0004 4 
0.008 7 0.0006 - 
“> 0.0067 <a. 0.0008 + 
—0.0010+ 
0.0044 
—0.0012+ 
nl 
9.004 —0.0014 + 
0.000 Porn 0.0016 
0 50 100 150 200 250 300 350 400 0 50 100 150 200 250 300 350 400 
t — [s] t — [s] 
0.010? 
0.008 + 
0.006 
0.0044 
0.002 - 
0.000 Porno 


0 50 100 150 200 250 300 350 400 0 50 100 150 200 250 300 350 400 
t — [s] t — [s] 


Figura 59: Modo Espiral 
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Figura 61: Modo Dutch Roll 
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Finalmente se determinaron los amortiguamientos, los periodos y tiempos de reducción a la 
mitad, los mismos se detallan en la Tabla[24] 


Tabla 24: Períodos, amortiguamientos , tiempos de reducción a la mitad y número de ciclos para 
reducción a la mitad. 


Modo Período Tiempo de reducción a la mitad Amortiguamiento Frecuencia 
[s] [s] ll [rad/s] 
Modo espiral - 49.4721 - - 
Modo de Rolido Puro - 0.2394 - - 
Modo Dutch Roll 2.5220 1.8032 0.1525 2.4914 


4.3.5. Modos trasversales de orden reducido 


= Modo espiral: Una buena aproximación para estimar la raíz del modo espiral es la siguiente: 


DA+E=0 
A, = =E[D 


E = 8 [(L,N; — LiN¡) 008 80 + (L,N; — LiN;) sin 60] 
D = =3 (£; cos 00 + N, sin 00) + uo (LN; — LN) 
= Modo de rolido puro: Una buena aproximación para obtener la raíz del modo de rolido 
puro es la propuesta por McRuer en 1973. Con este modelo se pueden obtener dos raíces 


que corresponden a aproximaciones para el modo de rolido y para el modo espiral. La apro- 
ximación es la siguiente: 


CÉ+DA+E=0 


C = UN, 
D= (LN; = Lo N;)= 8L, 
E= LN, = L.N,) 
= Modo de Dutch Roll: Para el modo de Dutch Roll, se pueda aproximar al valor propio 


mediante la siguiente expresión: 


B+BA+C=0 
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B= (Y, +N/) 
C= (YN” + UN ;) 


Resolviendo los modelos de orden reducido se obtuvieron los resultados presentados en la 
Tablas [25] y [26] Los valores propios obtenidos resultaron: 


Modo Espiral 1, = —-0,014049 


Modo Rolido Puro 1, = -3,106657 
Modo Balanceo del Holandés 2134 = -0,491465 + 2,408434¡ 


Tabla 25: Modelo de orden reducido para los modos espiral y de rolido puro 


Modo Espiral 
Parámetro Ecuaciones generales Modelo de orden reducido 
Amortiguamiento 0.0140 0.0141 
Tiempo 1/2 [s] 49.4721 49,3383 
Modo de Rolido Puro 
Parámetro Ecuaciones generales Modelo de orden reducido 
Amortiguamiento 2.8957 3.1067 
Tiempo 1/2 [s] 0.2394 0.2231 


Tabla 26: Modelo de orden reducido para el modo del Dutch roll 


Parámetro Ecuaciones generales Modelo de orden reducido 
Frecuencia amortiguada [rad/s] 2,4914 2,4084 
Período [s] 2,5220 2,6088 
Relación de amortiguamiento 0,1525 0,1999 
Tiempo 172 [s] 1,8032 1,4104 


4.3.6. Discusión 


Los resultados obtenidos indican que el comportamiento dinámico del avión es estable en 
relación con los modos principales, ya que en todo los casos, la parte real de los valores propios 
es negativa y como se aprecia en las gráficas existe una convergencia evidente alcanzada dentro de 
los plazos de tiempo esperados. Con la diferencia de que en el caso del modo Dutch Roll se tiene 
una solución convergente amortiguada, debido a sus raíces complejas conjugadas. 
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Cabe señalar que el Modo de Rolido Puro es el que tiene menor tiempo de convergencia ya 
que, como se esperaba, tiene el mayor módulo negativo y por lo tanto es altamente convergente. 

Analizando las Tablas[25]y[26] se puede concluir que los resultados derivados de las ecuaciones 
generales exhiben una notoria similitud numérica con los obtenidos mediante el modelo reducido. 
Esta semejanza abarca tanto los valores de amortiguamiento como los de tiempo reducido a la 
mitad, en el caso de los Modos Espiral y Rolido puro; y de la misma forma se extiende al Modo 
Dutch Roll en cuanto a los valores de frecuencia y periodo. Este hallazgo respalda la confiabilidad 
de los resultados, consolidando la validez y coherencia de ambos enfoques. 

De manera análoga al análisis realizado para el caso longitudinal, se examinó el caso trans- 
versal. El modo de rolido puro exhibe un amortiguamiento significativo. Al estudiar el modelo de 
orden reducido, se concluyó que a medida que el amortiguamiento de rolido C ,, aumenta en valor 
absoluto, el modo experimenta un incremento en el amortiguamiento. En nuestro caso, este valor 
puede ser considerable debido a que, en la estimación realizada mediante el método DATCOM, la 
pendiente de los perfiles alares del ala podría tener valores elevados, influyendo de manera signi- 
ficativa en los cálculos. En cuanto al Dutch Roll, la frecuencia del modo es alta, y en su modelo de 
orden reducido se observó que a medida que N, (relacionado con Cy,) aumenta, también lo hace 
la frecuencia del modo. 
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5. Conclusiones generales 


A lo largo de este estudio se ha llevado a cabo un análisis exhaustivo de la aeronave Short 360 
abordando distintos aspectos. 

En la primera etapa del estudio se analizó la performance de la aeronave, se determinaron 
distintos aspectos referido a esto y entre ellos se encuentran la determinación de la envolvente de 
vuelo de la aeronave, la determinación de los parámetros de ascenso como la tasa de ascenso 
RC y el gradiente de ascenso CGR, el análisis del comportamiento de la aeronave en vuelo en 
maniobra como sus velocidades de giro máximas y radios de giro mínimos, el estudio del alcance 
y autonomía y finalmente la determinación de las distancias de despegue y aterrizaje. Junto a 
esto fue posible verificar el cumplimiento de las normativas de aeronavegabilidad vigentes como 
también la comparación de los resultados obtenidos con lo provisto por el fabricante. 

En la segunda etapa, nos sumergimos en un análisis detallado de la Estabilidad Estática y Di- 
námica de la aeronave. En el ámbito de la Estabilidad Estática, abordamos tanto la estabilidad 
longitudinal como la lateral/direccional. En la primera, identificamos las posiciones extremas del 
C.G., trazamos la curva de sustentación en vuelo trimado y exploramos las relaciones entre el 
ángulo del elevador necesario para mantener el vuelo trimado y el ángulo de elevador requerido 
por g. Para la Estabilidad Estática Direccional y Lateral, determinamos las derivadas de estabili- 
dad mas importantes y determinamos la potencia del timón, el viento cruzado máximo admisible 
en el despegue, el vuelo con un motor inoperativo, la potencia de los alerones y la velocidad de 
rolido estacionaria. Luego, al adentrarnos en la Estabilidad Dinámica, estimamos las derivadas de 
estabilidad tanto en el ámbito longitudinal como en el lateral/direccional. Obtuvimos con éxito 
los modos longitudinales y transversales, y se pudo verificar y confirmar el cumplimiento de los 
criterios de estabilidad tanto en el dominio estático como dinámico. Este análisis exhaustivo pro- 
porciona una comprensión integral de la capacidad del Short 360 para mantener su estabilidad en 
diversas condiciones operativas. 

Para concluir, el análisis detallado de la aeronave Short 360 (Ver Figura [62) realizado a lo 
largo del curso ha sido una experiencia enriquecedora. La aplicación práctica de los conceptos 
aprendidos en la asignatura de Mecánica del Vuelo nos brindó una comprensión más profunda 
de los diversos aspectos que impactan en la performance y la estabilidad de una aeronave. Este 
estudio no solo cumplió con los objetivos académicos, sino que también resaltó la importancia de 
abordar la evaluación de aeronaves desde múltiples perspectivas. 


Figura 62: Aeronave Short 360 
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6. Apéndice 


6.1. 


La curva polar de resistencia de la aeronave se determinó para la condición de crucero y para 
la condición de baja velocidad, la mismas fueron obtenidas utilizando la metodología propuesta 


por Torenbeek (Torenbeek||1982). En la Figura[63|se muestran una gráfica con las curvas polares 


para 3 condiciones de vuelo. 


Curva polar 


6.1.1. Curva polar de alta velocidad 


Los resultados obtenidos por la metodología anteriormente mencionada son los que se mues- 
tran en la Tabla27] 


Tabla 27: Curva polar de alta velocidad 


Componente A B D 
Resistencia por sustentacion 
de los componentes 
Resistencia Ala sin alabeo 0.0000000  0.0000000  0.0269491 
Resistencia ala a sustentacion nula debida al alabeo 0.0000000  0.0000000  0.0000000 
Resistencia por sustentacion del fuselaje 0.0006767 -0.0024257  0.0021739 
Resistencia por sustentacion del empenaje 0.0007368  -0.0013129 0.0005848 
Resistencia de perfil 
Resistencia del ala 0.0057728  -0.0042194  0.0052743 
Resistencia del fuselaje - Upsweep 0.0099047  -0.0068000  0.0000970 
Barquillas 0.0029641  0.0000000  0.0000000 
Resistencia del empenaje horizontal 0.0002805  -0.0005000  0.0002230 
Resistencia del empenaje vertical 0.0026281  0.0000000  0.0000000 
Correcciones por interferencia 
Interferencia ala-fuselaje. Resistencia inducida por vortices 0.0004982  0.0000000  0.0000000 
Interferencia ala-fuselaje. Resistencia viscosa 0.0003350  -0.0003001  0.0000000 
Interferencia avion-empenaje 0.0000000 0.0003106  -0.0002767 
Barquillas 0.0000000  0.0069520  0.0000000 
Protuberancias, imperfecciones superficiales, etc 

Carenado tren principal 0.0007336 0.0000000  0.0000000 
Montantes ala 0.0001812 0.0000000  0.0000000 
Ala 0.0003464  -0.0002532  0.0003165 
Fuselaje + Empenaje 0.0000474  -0.0001698  0.0001522 
Superficies de control 0.0009315  -0.0001249  0.0000556 
Instalaciones 0.0000474  -0.0001698  0.0001522 
TOTAL 0.0260844  -0.0090132  0.0357017 
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6.1.2. Curva polar de baja velocidad 


Los resultados obtenidos por la metodología anteriormente mencionada son los que se mues- 


tran en la Tabla [28] 


Tabla 28: Curvas polares de baja velocidad (para despegue y aterrizaje) 


Curva polar de despegue con 0; = 15? 


Componente A B D 
Curva polar de alta velocidad 0.026084  -0.009013 0.035702 
Resistencia de perfil 0.003275  -0.003207 0.000000 
Resistencia vorticosa 0.011199  -0.003151  0.000000 
Resistencia de trimado 0.001368 0.000033 —0.000000 
Resistencia por tren de aterrizaje 0.009209 -0.003357 0.000306 
TOTAL 0.051135  -0.018695 0.036008 
Curva polar de aterrizaje con 6; = 309 
Componente A B D 
Curva polar de alta velocidad 0.02608  -0.00901  0.03570 
Resistencia de perfil 0.00864 — -0.00552  0.00000 
Resistencia vorticosa 0.01929 — -0.00543 — 0.00000 
Resistencia de trimado 0.00270  0.00007  0.00000 
Resistencia por tren de aterrizaje  0.00229  -0.00167  0.00031 
TOTAL 0.05900  -0.02157  0.03601 
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0.164. —— Polar condición crucero ' 
—— Polar df = 15” 


—— Polar 0; = 30” 
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Figura 63: Curva polar de resistencia para tres condiciones de vuelo 
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6.2. Curva de potencia de los motores 


La aeronave está equipada con motores Pratt 8 Whitney PT6-63R de la serie PT6-60. Tanto en 
el manual del fabricante de la aeronave como en los manuales del fabricante del motor, se pueden 
encontrar curvas, como se muestra en la Figura[64](Ver Ref. 2007)). 

Estas curvas indican que, para motores de la serie PT6-60, el torque en [/bs— ft] proporcionado 
por el motor es constante a una presión atmosférica específica, expresada en pulgadas de mercurio 
[in — Hg]. Conforme aumenta la altitud y disminuye la presión atmosférica, el motor es capaz de 
entregar un porcentaje reducido de su torque máximo, como se espera de acuedo a lo visto en la 
teoría. 

Utilizando la información proporcionada por el fabricante del motor acerca del valor máximo 
de torque que el motor puede generar, hemos podido derivar una expresión que relaciona el torque 
del motor con la altitud (Ver Expresión 79). 


T = 0,0897 - h + 3625 [Ib - feet] (79) 


Para calcular la potencia en el eje del motor, necesitamos tener en cuenta las RPM especificadas 
por el fabricante de la aeronave, que son 1300 RPM. Dado que la aeronave tiene dos motores, 
debemos considerar este factor. Además, al realizar las conversiones adecuadas para pasar de [/bs— 
ft] a [Vm] y convertir las RPM a [rad/seg], podemos obtener la potencia en el eje del motor como 
se muestra en la Expresión|80; 


2 
Peje =2-T E -1,3558 - 1300) [Watts] (80) 


en donde 


z, Cantidad de motores 

factor de conversión para pasar de RPM a [rad/seg] 

1.3558 factor de conversión para pasar de [lbs — ft] a [Nm] 
1300 Revoluciones por minuto del motor [RPM] 
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Figura 64: Curvas pertinentes a los motores PT'6-60s 
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6.3. Curva de eficiencia de la hélice 


Dado que no se disponía de curvas proporcionadas por el fabricante, se optó por emplear el 
método numérico propuesto por U.P Solies en su trabajo titulado "Numerical Method for Estima- 
tion of Propeller Efficiencies"(Solies|[1994). Los resultados de este método se obtienen mediante 
la resolución de un proceso iterativo regido por las siguientes expresiones: 


n= Ho Mi 


T = (97 + P/V) 


Eo data 

2 2 4 2pA 

Ñ V 

V+ (2%) 

Este método proporciona una estimación fundamental de la eficiencia de la hélice cuando no 


se dispone de información directa del fabricante. Los resultados obtenidos se detallan en la Tabla 


Mi 


Tabla 29: Resultados obtenidos. 


Vo on 60.00 0.76 
10.00 0.53 65.00 0.77 
15.00 0.58 70.00 0.78 
20.00 0.62 75.00 0.79 
25.00 0.65 80.00 0.80 
30.00 0.67 85.00 0.80 
35.00 0.69 90.00 0.81 
40.00 0.71 100.00 0.82 
45.00 0.73 110.00 0.83 
50.00 0.74 120.00 0.84 
55.00 0.75 150.00 0.86 


Luego, se realizó un ajuste numérico a los puntos obtenidos utilizando las herramientas pro- 


porcionadas por la librería Scipy de Python (consultar scipy.optimize.curve_fit). La expresión de 


ajuste se presenta en la Ec. Los resultados y la curva ajustada se muestran en la Figura 


1 = -2,7837E-09-V*+1,0935E-06-V*+-1,6159£E-04-V?+1,1873E-02-V+4,3272E-01 (81) 
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MW Resultados Solies 


=— 1), Ajuste con polinomio de cuarto grado 


0 20 40 60 SO 100 120 140 160 
V [m/s] 


Figura 65: Representación de la eficiencia de la hélice. 
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6.4. Consumo específico en crucero 


A partir de la Figura[64]fue posible relevar la variación del flujo de combustible con la presión 
atmosférica (la cual se encuentra dada en pulgadas de mercurio). Utilizando las expresiones de la 
atmósfera estándar, teniendo en cuenta los dos motores y realizando las conversiones de unidades 
correspondientes se pudo ajustar una expresión de flujo de combustible [V/h] vs altitud [m]. Dicha 
expresión resulta: 


Fp = -0,6424 - h + 5621,40 


Obtenida una expresión para la variación del flujo de combustible, se obtiene el consumo es- 
pecífico en crucero haciendo uso de la Expresión |S0|como sigue: 


Fr 
CA = 
Pp 

Peje 


Finalmente el consumo específico para el crucero resultó: 


Cerúcero E 1 > 17E-06 ES 
Ss 
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6.5. Velocidades de pérdida y coef. sustentación máx. 


6.5.1. Condición de despegue 


Para la condición de despegue, se determinó la velocidad de pérdida considerando el peso 
máximo al despegue (12000 kg) y una deflexión de flaps de 15”. Con esta información, se pudo 
determinar la velocidad de pérdida consultando el manual del fabricante de la aeronave, como se 
muestra en la Figura|66] 


MINIMUM SPEED IN STALL - POWER OFF 


D.A.5B43.P2155 


STALL SPEED - KNOTS C.A.S. 


17 18 19 20 21 22 23 24 25 26 27 
AIRCRAFT WEIGHT - LB/1000 


Figura 66: Velocidad de pérdida en condición de despegue 


Dado que la velocidad de pérdida obtenida del gráfico en la Figura [66] corresponde a la velo- 
cidad calibrada (CAS), se procede a calcular la velocidad verdadera (TAS). Es importante señalar 
que se asume un factor f. igual a uno, ya que no se tienen en cuenta los efectos de compresibilidad. 


EAS = f.:CAS =TAS - 


m 
TAS stall,TO = Vstall.TO = 44,24 E 


Luego por otro lado, teniendo ya esta velocidad se determinó el C,, para la condición de 


despegue, el mismo resulta: 
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- MM 224333 


2 
PS V calro 


Emax 


6.5.2. Condición de aterrizaje 


Para la condición de aterrizaje se procedió de forma análoga a lo detallado anteriormente, para 
la condición de aterrizaje se tomó el máximo peso al aterrizaje (11840 kg) y una deflexión de flaps 
de 30”. Con esta información, se pudo determinar la velocidad de pérdida consultando el manual 
del fabricante de la aeronave, como se muestra en la Figura 


MINIMUM SPEED IN STALL - POWER OFF 


D.A.SB43.P2155 


STALL SPEED - KNOTS C.A.S. 


AIRCRAFT WEIGHT - LB/1000 


Figura 67: Velocidad de pérdida en condición de aterrizaje 
Luego resultaron: 


m 
E A E 


2W 
Cha = === — = 2,7065 
ps v? 


stall,LD 
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6.6. Pendiente ala - fuselaje 


Para la determinación de la pendiente de sustentación ala - fuselaje (C,, ),, se siguió la meto- 
dología propuesta por la cátedra de Aerodinámica I (Cátedra de Aerodinámica I| 2021). La misma 
se encuentra dada por la Expresión[82] Basicamente consiste en multiplicar a la pendiente del ala 
por un factor de corrección K;. 


(Cr wf = Mw» = KAC,) yw (82) 


donde el factor K, se calcula de la siguiente forma: 


by So T bi 
K¡=1+2,15—-=+ PE 
E DS» ACL)w Sy 


en donde: 


(CL,)w = 5,7597 Pendiente del ala [1/rad] 
S, =42,73 Superficie alar [m?] 
S.=36,79 Superficie alar expuesta [m?] 
b= 22,80 Envergadura [m] 

by =2,30 Ancho del fuselaje [m] 


El valor de K, = 1,0815; por lo que la pendiente del ala - fuselaje resultó: 


Awp» = 6,2290 [1 /rad] 
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6.7. Centro aerodinámico ala - fuselaje 


La posición del centro aerodinámico de la combinación ala-fuselaje se puede obtener aplicando 


la Expresión [83|(Cátedra de Aerodinámica 1)/2021). 


ac ac A ac A ac 
(Je) = (ze) y Alte , Afetoe (83) 
C Iwf C /w C Cc 
en donde: 
AfiXac 18 Dyhyly, 
Cc - (Cr. wp Se 


es una corrección experimental de la contribución de las secciones aguas arriba del ala y 


AfoXac 0,273 bfcg(b—by) 
¿E (1+2)22(b+2,15bj) 
Tiene en cuenta la pérdida de sustentación del ala dentro del fuselaje y es derivada del método 


DATCOM. 
en donde: 


tan Aja 


(Cr.)wf = 5,7597 Pendiente del ala [1/rad] 
by = 2,30 Ancho del fuselaje [m] 
hy = 2,82 Alto del fuselaje [m] 


ls, = 8,60 Distancia nariz del fuselaje al B.A. [m] 
c=1,99 Cuerda media aerodinámica [»m] 
Cg = 1,85 Cuerda media geométrica [m] 


s$=142.73 Superficie alar [m?] 
b= 22,80 Envergadura [mm] 
Ay /4c = 0,30 Ángulo de flecha al cuarto de cuerda [DEG] 


Resulta: 


Xac 


(E) = 0,0611 (84) 
C Iwf 
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6.8. Coeficiente de momento libre ala - fuselaje 


Para obtener el valor del coeficiente de momento libre de la configuración ala-fuselaje se siguió 
la metodología propuesta por la cátedra de Aerodinámica I ((Cátedra de Aerodinámica 1//2021)), 
representada por la Expresión [85] 


(A E (Ora. )., + AFC mie (85) 


Donde (C,,.),, corresponde al momento libre del ala y la corección AFC, se obtiene a partir 
de la teoría de Munk, con la Expresión [86] 


2.10s Tb gh pl y Cr 
AfCm. ==1,8|1 — 86 
a 09 


Teniendo como datos: 


(Como )w = 0,0805 Coeficiente de momento libre del ala [Adim] 


bs = 2,30 Ancho del fuselaje [m] 

hy = 2,82 Alto del fuselaje [m] 

If = 21,2 Largo del fuselaje [m)] 

c=1,99 Cuerda media aerodinámica [m] 
Ss =4078 Superficie alar [m?] 


(Cr.)wf = 5,7597 Pendiente del ala [1/rad] 
Ch, = 0,4866 Coeficiente de sustentación libre [Adim] 


Obteniendo así: 
(Crmac) ob => 0, 2106 
Análogo al procedimiento anterior, pero teniendo en consideración la condición de despegue 
del aeronave, se obtiene: 


(Cors) = 0, 4534 
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6.9. Pendiente del empenaje horizontal y el empenaje vertical 


6.9.1. Empenaje Horizontal 


Para determinar la pendiente de sustentación del empenaje horizontal se empleó la formula 
de Polhamus (Ver Expresión [87), en dicha expresión se requiere el alargamiento Az del empenaje 
horizontal, el parámetro de Prandtl Glauert BG, la flecha a la mitad de cuerda A, y la relación 
entre la pendiente de sustentación del perfil y la téorica de perfiles delgados k. 


2TAR 


2 2 
DE ArBPG a [o +4 
k P 


Ch (87) 


ah 


PG 


en este caso Ag = 5,6918, la pendiente del perfil NACA 0012 c,, = 5,7869 -[1/radl, Bpú = 0,9525 
y A¿¡2 = 07. La pendiente del empenaje horizontal resulta: 


Ch,, = Y = 4,3523 [1/rad] 


6.9.2. Empenaje Vertical 


Para la determinación del empenaje vertical también se empleó la formula de Polhamus (Ver 
Expresión|87) aunque previamente para aplicarla se determinó el alargamiento efectivo del empe- 


naje vertical, Az, ,,,, este se determina a partir de la Expresión |88] 


A Ms 
Also = 2) [ + Ky [Green a 1) =3,1770 (88) 


Ar, Ry.B) 


en donde los parámetros necesarios para aplicar la Expresión [88|se obtuvieron de las gráficas de 
las Figura Luego aplicando la formula de Polhamus con la pendiente del perfil NACA 0016 
Cc, = 5,5290 - [1/rad], BpúG = 0,9525 y Ay, = 22”. La pendiente del empenaje horizontal resulta: 


Ch, , = 4y = 4,2288 [1/rad] 


e, 
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Figura 68: Parámetros para calcular Ar, ,,, 
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6.10. Gradiente de Downwash 


Para determinar la expresión general del downwash promedio en el empenaje horizontal del 
avión asignado utilizando el método ESDU 80020 (ESDU|/2012). El método es aplicable a confi- 
guraciones dentro de los límites siguientes: 


Alargamiento 2<A<12 
Ahusamiento 0<2<l 
Flecha 0 < Aya < 52" 
Ancho del fuselaje 0,1 <w/b < 0,2 
Posición € = 2 0,5 <ér <2 
Posición 77 = 55 0,2 < nr < 0,8 
Posición ¿7 = 53 0,2 < £y <0,6 
Número de Reynolds basado en¿ 0,5 x 10% < Res < 15 x 10% 
Número de Mach Los que permitan las reglas de similitud subsónica 


6.10.1. Proceso de cálculo 


El procedimiento se encuentra sintetizado en el siguiente listado: 


1. 


Con el alargamiento A y el producto de A - tan(A;y2) en la Figura 13 del report (Ver Figuras 
del report) se obtiene Hp... 


. Con la Tan(A;4) y la inversa de 1/£7 en la Figura 14 (Ver Figura del report) se obtiene Hp,. 


Empleando la Expresión |89|se determina Hp. 


Hp = Hp, — (1 Hp.) (89) 


. Con el alargamiento A y Hp, en la Figura 15 (Ver Figura del report) se obtiene AH» 


. Con el ahusamiento 2, el alargamiento A, 77 y A : Tan(A;,) en las Figuras 16, 17,18 y 


19 (Ver Figuras del report), se obtiene F. Mediante la Expresión [90] se puede obtener el 
gradiente promedio del downwash en el empenaje cuando éste está inmerso en la lámina 
vorticosa (¿7 = 0). 


deYy 2a 
—| = —-—(H AHp)F 90 
1) Al p + AHp) (90) 


. Con el alargamiento A en Figura 20 (Ver Figura del report) se puede obtener B. Con la 


Expresión |91|se puede obtener la variación del downwash en el plano de la estela. 


dé dev 
(de), Bla) pe 
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2(L) 
Oy MXda 
<F 


de 
6. Con (o en la Figura 21 se obtiene 
A 


UN (£) 
O£r Xda 


Para Mach = 0 se puede contemplar el efecto de la distancia entre ambos mediante la 


Relación [92] 
de de O (de 
A E 
00 da), 067 | da 


7. Con el valor del Mach en crucero, en la Figura 23 (Ver Figura del report) se obtiene Ky y es 
posible corregir la expresión[92|con la expresión [93] 


de de 
e), (7), o E 


8. Con la diferencia enter ¿7 — ¿g, en la Figura 22 (Ver Figura del report) se obtiene €;,. Final- 
mente con la expresión [94]se obtiene la expresión para e. 


por lo que se determina exactamente el valor de 


; para colas bajas se debe multiplicar este factor por 1,3. 


(92) 


E € 
€ = El, + [dj — Q11p) (94) 
Q 
Para la aeronave luego de aplicar la metodología esto resultó: 


€ = 0,0597 + 0,3640 a 
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6.11. Momentos de Inercia 


Los momentos de inercia de la aeronave fueron estimados siguiendo la metodología propuesta 


por la bibliografía recomendada por la cátedra (Chawlaj|1952). 


Los momentos de inercia fueron estimados teniendo en cuenta los radios de giro adimensio- 
nales de referencia para el tipo de aeronave que estamos tratando, para este caso, la misma es una 
aeronave de transporte. Los valores tomados como referencia fueron los dados por la Tabla [30] 


Tabla 30: Radios de giro adimensionales para la aeronave 


Radios de giro adimensionales 
Tipo de avión R, R, R, 
Carga - Transporte 0,2900 0,3300 0,4220 


Los radios de giro adimensionales R,, R, y R; se encuentran definidos por la Expresión[95] 


Kx 
R,= 
b/2 
k 
R) == 95 
k 
R,= 
¡e 


Luego, a partir de ellos, se pueden determinar los radios de giro dimensionales teniendo en 
cuenta los valores de la envergadura b = 22,80 [m], el largo de la aeronave L = 21,58 [m] y 
el parámetro e = (1/21(b + L) = 22,19 [m]. Los valores de dichos radios de giro dimensionales 
resultan: 


k, = 3,3060 [m] 
k, = 3,5607 [m] 
k. = 4,6821 [m] 


Finalmente a partir de los resultados obtenidos, se pueden estimar los momentos de inercia 
siguiendo la Expresión[96| en donde W/g = 9587 [kg] (masa en condición crucero). 


ie (96) 
8 


Los momentos de inercia de la aeronave resultan: 
L, = 104785 [Kgm”] 


1, = 121553 [Kgm/] 
L = 210172 [Kgm?] 
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